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Wofi:  Completed 

A.  Implicit  Clouds: 

We  have  made  considerable  progress  in  the  area  of  implicit  treatment  of  clouds.  This 
iminovement  resulted  partly  from  involving  physical  initialization  prior  to  the  global  model 
forecast  experiments.  A  review  of  this  woik  appears  in  Appendix  1  of  this  proposal. 

a)  Physical  Initialization: 

When  we  examine  the  frequencies  of  a  complete  atmospheric  nonhydrostadc  model  we 
find  that  acousdc  waves,  gravity  and  slow  moving  atmospheric  Rossby  waves  are  the  principal 
modes.  Errors  in  the  inidal  data  analysis,  especially  of  the  mass  and  motion  field  result  in  the 
growth  of  the  higher  frequency  motions.  We  call  these  computational  modes.  The  use  of 
hydrostatic  approximation  in  the  vertical  equation  of  motion  is  adequate  to  suppress  the 
acoustic  waves.  The  suppression  of  high  frequency  computational  gravity  modes  is  usually 
accomplished  using  a  normal  mode  iititialization. 

Given  an  initial  data  set  which  is  subjected  to  an  analysis  scheme  (such  as  optimal 
interpolation  or  variational  statistical  methods),  an  assimilation  within  the  model's 
infrastructure  and  to  the  normal  mode  iititialization  one  would  expect  to  have  described  a 
reasonable  initial  state  for  the  model's  forecast  Such  is  not  the  case  if  there  are  major  data 
gaps  as  is  the  case  especially  over  the  tropical  latitudes  and  the  southern  hemisphere.  There 
exist  some  alternate  data  sets  such  as  measures  of  rain  rates  from  remote  sensing  from 
satellites.  Rainfall  rate  is  not  a  basic  dependent  variable  of  a  forecast  model.  It  is  an  internal 
product  To  incorporate  rainrates  into  the  data  analysis  one  needs  to  assimilate  a  precipitation 
process  consistent  with  the  rainrates  as  input  That  can  be  done  using  inverse  algorithms  that 
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^nerate  proxy  variables  for  the  assimilation.  This  method  is  called  physical  initialization, 
Krishnamurti  et  al.  (1991a). 

The  humidity  variable  is  poorly  described  by  the  available  observations,  especially  over 
the  tropics.  Physical  initialization  within  a  data  assimilation  restructures  the  vertical 
distribution  of  humidity  consistent  with  the  observed  rainfall  rates.  This  enables  the  model  to 
return  a  rain  rate  nearly  the  same  as  the  ob^rved  during  the  physical  initialization  within  the 
data  assimilation.  During  the  assimilation  the  heating  from  the  condensation  process,  the 
divergence  distribution  and  the  surface  pressure  fields  adjust  to  the  rain  rates.  This  procedure 
has  been  found  to  exhibit  a  considerable  skill  in  the  rainfall  forecasts  over  very  short  range 
time  scales.  The  noteworthy  improvements  are  in  the  nowcasting  (day  0)  and  the  one  day 
forecast  of  rainfall  over  the  entire  tropics.  Figure  (la)  illustrates  a  24  hour  tropical  rainfall 
distribution  as  measured  by  satellite  (microwave  radiometers  and  outgoing  longwave  radiation) 
and  raingauge  networic,  and  figure  (lb)  illustrates  the  model  initialization  of  the  24  hour  rain. 
As  is  evident  here  a  high  resolution  global  model  (Triangular  truncation  170  waves)  is  able  to 
initialize  such  rainfall  very  closely.  These  rainfall  observations  supplement  the  other 
conventional  weather  observations  via  the  inverse  algorithms.  The  day  0  and  the  one  day 
forecast  skill  of  rainfall  is  illustrated  in  figure  (2a,b).  Here  we  show  correlation  of  observed 
and  initialized  (or  predicted)  rainfall  for  data  over  transform  grid  squares.  Also  shown  in  these 
illustrations  are  some  correlations  for  operational  weather  centers  and  for  our  control 
experiments  where  such  a  physical  initialization  was  not  executed.  Overall  it  is  apparent  that 
the  rain  rate  initialization  provides  a  very  useful  data  source  for  the  tropical  latitudes.  In  the 
appendix  we  show  a  strong  positive  impact  on  cloud  forecasts  from  physical  initialization. 

b)  Definition  of  Implicit  Clouds: 

Clouds  are  classified  as  high,  middle,  and  low.  High  clouds  are  defined  to  exist 
between  400  and  100  mb,  middle  clouds  between  700  and  400  mb,  and  low  clouds  between 
900  and  700  mb.  Clouds  are  not  defined  for  sigma  levels  below  low  cloud  bases  or  above 
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high  cloud  tops. 

The  amount  of  a  cloud  type  in  an  individual  layer  is  assumed  to  be  a  function  of  the 
relative  humidity,  which  is  dependent  on  the  moisture,  temperature,  and  pressure.  The  clouds 

are  assumed  to  be  present  when  the  mean  relative  humidity  (RH)  in  a  layer  exceeds  the 
threshold  value  (RHc).  The  cloud  amount  N  for  each  type  of  cloud  is  defined  by  the  relation. 


N 


(h,  m  or  l) 


'  RH  -  RHc 
1  -  RHc 
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RH>RHc 


where  RHc  has  been  set  to  0.25,  0.40,  and  0.50  for  high,  medium,  and  low  clouds, 

respectively.  RH  is  the  pressure  weighted  mean  relative  humidity  in  a  layer.  The  maximum 

possible  value  of  N  is  1.  The  cloud  amount  N  is  set  to  0,  when  RH  is  less  than  RHc.  Total 
cloud  amount  (Np  is  determined  by  using  random  overlap  of  the  cloud  layers  in  a  column  as 
given  by  =  1  -  (l-Nj^)(I-Nj^)(I-Njj).  Cloud  covers  (CQ,)  is  partitioned  into  eight  sky 
condition  categories  as  defined  in  Kiishnamurti  et  al.  (1990); 


CCi  =  (l-NJ(l-Nj^)(l.Njj) 
CC2  =  Nj^(1-N^)(1-Njj) 

CC3=N^N„(1.N„) 

CC5  =  N^(1-N^)(1-N^) 

CCs  =  NmN^(1-N^) 

cx:7  =  n^(1-n^)(i-n^,) 

CC.  =  N,N„(1-NJ 


clear  sky 
low  clouds  only 
low  and  middle  clouds 
low,  middle,  and  high  clouds 
middle  clouds  only 
middle  and  high  clouds 
high  clouds  only 
low  and  high  clouds 


This  sky  condition  is  used  in  the  computation  of  the  solar  radiation  absorbed  or 
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scattered  by  clouds. 


c)  Treatment  of  Cirrus  Clouds: 

High  level  clouds  (e.g.  cirrus)  play  an  important  role  in  the  atmospheric  radiation 
budget  They  allow  bodi  the  short  and  longwave  radiation  to  be  transmitted  through  them 
partially,  and  therefore  are  treated  as  non-black  in  the  radiative  transfer  model.  Their 
characteristics  are  based  on  the  optical  depth,  the  single  scattering  albedo  and  the  symmetry 
scattering  factor  in  the  radiation  scheme.  The  optical  depth  of  cirrus  is  simply  assumed  to  be  a 
function  of  temperature.  The  contents  of  cirrus  categorize  ice  and  supercooled  water  droplets 
with  temperature  range  between  190.65  K  and  263.15  K,  supercooled  water  droplets  between 
263.15  K  and  273.15  K.  The  single  scattering  albedo  of  cirrus  is  then  determined  based  on  the 
ratio  of  the  optical  depth  of  cirrus  to  the  total  optical  depth  of  the  atmosphere.  The  symmetry 
scattering  factor  is  fixed  as  constant  for  clouds. 

No  cloud  microphysics  based  on  phase  change  of  water  is  considered  in  the  implicit 
cloud  scheme. 

d)  Treatment  of  Tropical  Clouds: 

In  general,  cumulus  clouds  have  horizontal  scales  of  the  order  of  a  few  kilometers,  and 
are  not  explicitly  resolved  by  large  scale  numerical  weather  prediction  models.  Cumulus 
parameterization  was  needed  to  address  this  issue.  Our  cumulus  parameterization  for  deep 
cloud  in  a  tropical  region  is  based  on  Kuo  (1974).  The  heating  and  moistening  of  the 
atmosphere  and  the  rainfall  rates  of  this  procedure  are  simple  functions  of  the  large  scale 
convergence.  A  further  refinement  of  this  scheme  was  provided  by  Krishnamurti  and  Bedi 
(1988a).  Here  statistically  based  improvements  on  the  Kuo  scheme  were  made  to  provide 
global  applicability  and  minimization  of  initial  errors  of  the  heating,  moistening  and  rainfall 
rates. 

For  convective  clouds,  a  cloud  fraction,  CC,  is  based  on  the  rainfall  rates  of  a  Kuo's 
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scheme.  The  past  three  hours  of  rainfall  total  are  used  to  determine  the  cloud  coverage 
firacdon  CC.  Table  6  describes  these  relationships. 

Table  5:  Relationship  Between  the  Precipitation  Rates  and  Convective  Cloud 

Cover  in  NMC  Model 

P(mm/day)  .14  .31  .70  1.6  3.4  7.7  17.  38.  85. 

CC  (fraction)  0  .1  .2  .3  .4  .5  .6  .7  .8 

For  rainfall  rates  less  than  1.6  mm/day,  25%  of  the  cloud  coverage  fraction  for 
convective  clouds  are  identified  as  high  clouds  and  75%  as  the  low  cumulus  bed.  However,  if 
the  rainfall  exceeds  1.6  mm/day,  then  the  cloud  fraction  for  cirrus  is  given  by  a  relation, 

C.  =2*(CC-.3) 


This  is  supposed  to  account  for  the  cumulonimbus  anvils. 

e)  Treatment  of  Radiative  Transfer: 

In  most  models  the  physics  algorithms  include  rather  sophisticated  parameterizations 
for  the  radiative  transfer  processes.  The  total  problem  includes  the  parameterization  of  short 
and  long  wave  irradianccs,  the  effects  of  clouds,  role  of  atmospheric  constituents  such  as  water 
vapor,  carbon  dioxide  and  ozone.  The  definition  of  clouds  and  cloud  radiative  interactions  are 
important  areas  in  modelling.  The  surface  energy  balance  is  also  strongly  radiatively 
controlled  since  it  effects  issues  such  as  the  diurnal  changes  and  land  surface  processes. 

Most  current  radiative  transfer  models  include  what  are  called  band  models.  The 
atmospheric  absorption  bands  for  the  solar  and  the  long  wave  irradiances  are  sorted  out  over 
selected  representative  absorption  bands.  This  is  a  time  saving  technique  and  it  is  possible  to 
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replicate  results  (for  clear  sky  conditions)  as  one  would  obtain  from  line  by  line  integrations 
across  the  entire  atmo^heric  spectrum. 

For  clear  skies  we  use  the  so-called  k-distribution  approximation  following  Stephens 
(1^)  for  the  absorption  of  solar  radiation.  For  cloudy  sides  we  invoke  a  detailed  multiple 
scattering  computation  based  on  the  work  of  Davies  (1982).  Here,  apart  from  absorption,  we 
consido'  tile  reflection  and  transmission  characteristics  of  various  layers  of  the  atmosphere. 
Following  Joseph  and  Wiscombe  (1976),  we  use  the  so-called  Delta-Eddington  approximation 
which  is  also  called  a  two  stream  approximation.  The  ozone  absorption,  especially  in  the 
strato^here,  is  important  where  the  Rayleigh  scattering  is  significant  This  is  mostly  below 
tile  ozone  layer  of  the  stratosphere.  Consideration  is  given  to  the  weak  visible  bands  (Chapius 
bands)  and  the  ultraviolet  bands  (Hartley  and  Higgins  bands)  following  the  formulation  of 
Lacis  and  Hansen  (1974).  Surface  absorption  is  divided  into  two  bands,  one  where  the  water 
vapor  absorption  is  large  and  the  other  includes  the  rest  of  the  spectrum.  The  parameterization 
of  long  wave  irradiances  follows  Harshvardhan  and  Corsetti  (1984).  This  includes  band 
scaling  approximations  (the  band-centers  and  band-wings)  for  carbon  dioxide  and  water  vapor 
following  Chou  and  Peng  (1983)  and  Chou  (1984).  The  transmission  functions  for  ozone  are 
based  on  the  formulation  of  Rodgers  (1968).  The  molecular  line  absorption  and  e-type 
absorption  (continuum  absorption)  of  diffuse-transmittance  function  due  to  both  water  vapor 
and  CO2  ^  included  in  this  analysis. 

The  specification  of  clouds  in  many  models  is  based  on  certain  threshold  humidity 
criteria.  Three  types  of  clouds  are  considered:  low,  middle  and  high.  If  the  model  based 
predicted  humidity  exceeds  certain  thresholds  then  the  radiatively  active  clouds  are  assumed  to 
exist  in  that  vertical  layer  of  the  atmosphere.  Low  and  middle  clouds  are  regarded  as  black 
body  radiators;  the  high  level  clouds  (e.g.  cirrus),  however,  are  treated  as  partial  black  body 
radiators. 

An  example  of  cloud-radiative  interaction  is  presented  next  Tropical  oceans  exhibit  a 
very  large  coverage  of  shallow  stratocumulus  clouds.  The  atmosphere  over  the  tropical  ocean 


6 


exhibits  a  moist  instability  over  the  lower  troposphere.  The  cloud  base  over  this  region  lies 
roughly  1  km  above  the  ocean  and  the  cloud  top  is  usually  located  near  the  top  of  this  unstable 
region  \^ch  is  roughly  3  km  above  the  ocean.  The  maintenance  of  this  unstable  layer  is  a 
central  issue  in  tropical  meteorology.  Radiative  processes  and  surface  evaporation  contribute 
to  the  steady  restoration  of  this  instability  whereas  convective  clouds  contribute  towards  its 
erosion.  Tropical  weather  disturbances  feed  on  tliis  instability.  If  the  restoring  mechanisms 
are  not  properly  modeled,  the  air  entering  tropical  storms  does  not  provide  sufficient  instability 
for  cloud  growth  and  storm  forecasts  tend  to  degrade.  The  long  wave  radiative  cloud  top 
cooling  is  one  of  the  major  mechanisms  that  restores  conditional  instability.  Fig.  (3a,  b) 
illustrates  results  of  4  day  forecasts  of  a  storm  where  two  different  radiative  parameterization 
methods  were  used  in  the  global  model.  The  bottom  panel  shows  the  low  cloud  distributions 
for  the  above  band  model.  The  upper  panel  illustrates  the  results  from  a  classical  radiation 
algorithm  called  the  emissivity-absorptivity  method,  Chang  (1980).  The  latter  did  not  provide 
a  sufficient  cooling  rate  over  the  lower  troposphere;  the  temperatures  were  high,  the  relative 
humidity  as  a  consequence  was  too  low,  and  low  clouds  did  not  form  since  the  threshold  value 
of  relative  humidity  were  not  exceeded.  Lacking  sufficient  low  clouds  the  restoration  of 
conditional  instability  was  underestimated.  Hg.  (4)  illustrates  the  time  histories  of  the 
predicted  cooling  rates  and  humidity  for  these  two  respective  radiative  transfer  models.  In 
these  experiments  a  vigorous  storm  did  form  and  move  over  China  when  the  band  model  was 
used,  and  that  storm  weakened  and  disappeared  in  the  forecast  where  the 
emissivity-absorptivity  model  for  cloud  radiative  interaction  was  used,  Krishnamurti  et  al. 
(1991b). 

f)  Results  from  the  implicit  cloud  scheme: 

These  are  presented  in  Appendix  1. 

B.  Explicit  Clouds: 

We  have  made  some  progress  in  the  areas  of  the  explicit  treatment  of  clouds.  This  has 
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required  a  major  extension  of  the  global  model  from  the  inclusion  of  two  dependent  variables 
i.e.,  the  cloud  water  mixing  ratio  and  the  cloud  fractions.  This  work  is  still  ongoing  and  will 
form  die  major  area  of  research  for  the  new  proposal.  Cloud  microphysical  processes  are 
crudely  represented.  In  this  scheme  the  time  evolution  of  the  cloud  variables  are  determined 
by  the  sources  and  sinks  due  to  the  several  cloud  processes.  This  scheme  considers  the 
formation  of  clouds  in  connection  with  large-scale  ascent,  cumulus  convection,  and  radiative 
cooling.  Clouds  dissipate  through  adiabatic  and  diabatic  heating,  turbulent  mixing  of  cloud  air 
with  unsaturated  environmental  air,  and  depletion  of  cloud  water  by  precipitation.  Both  liquid 
and  ice  phases  are  included  in  the  cloud  model.  The  liquid  phase  processes  consist  of 
evaporation,  condensation,  autoconversion  and  precipitation.  The  ice  phase  processes  include 
only  gravitational  settling  of  cloud  ice.  Cloud  water  and  cloud  ice  are  presently  distinguished 
based  on  the  prescribed  temperature  but  are  treated  as  a  bulk  for  cloud  water  path  in  the 
radiation  scheme.  Cloud  processes  are  represented  using  a  bulk  water  parameterization 
technique,  where  the  liquid  phase  is  subdivided  into  cloud  water  and  precipitation  water. 
Cloud  water  is  assumed  to  form  when  the  relative  humidity  exceeds  a  specified  threshold 
value,  and  precipitation  processes  are  described  by  simple  empirical  parameterizations  fcr 
warm  and  cold  clouds. 

The  advantage  of  this  approach  is  that  it  is  physically  realistic  and  is  consistent  with 
the  rest  of  the  model,  while  directly  accounting  for  the  various  feedbacks  involving  clouds. 
However,  the  accuracy  of  cloud  parameterization  within  this  approach  depends  critically  on  the 
realistic  treatment  of  advective  transport  of  cloud  variables,  subgrid-scale  processes,  cloud 
microphysics,  and  cloud  optical  properties.  Note  that  there  still  exist  large  uncertainties  in  the 
representation  of  these  processes.  Details  of  these  are  presented  in  Appendix  1. 
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ABSTRACT 


The  impact  of  the  initial  data  on  the  cloud  forecast  by  the  Florida  State 
University  Global  Spectral  Model  (FSUGSM)  has  been  investigated.  Additionally  two 
different  cloud  schemes,  one  implicit,  the  other  explicit,  were  tested.  This  work  has 
shown  that  improving  the  information  content  of  the  initial  data  by  physical 
initialization  has  a  very  strong,  positive  impact  on  cloud  forecasts. 

The  physically-initialized  data  has  a  greater  impact  in  the  implicit  cloud  scheme. 
Model  spin-up  of  cloud  generation  is  considerably  reduced.  There  is  an  overall  better 
representation  of  hi^,  middle,  low  and  total  clouds  over  the  tropics  and  there  is  a 
discernible  improvement  in  the  prediction  of  clouds.  A  strong  negative  correlation 
between  cloud  shortwave  forcing  and  longwave  forcing  has  been  observed  in  the 
implicit  scheme  experiments.  This  result  matches  very  closely  with  the  Earth  Radiation 
Budget  Experiment  (ERBE)  observations.  The  explicit  scheme  docs  not  respond  in  a 
similar  fashion.  The  physically-initialized  data  has  only  a  minor  impact  on  forecasts 
made  using  the  explicit  scheme.  This  suggests  a  need  for  further  research  into  the 
physical  initialization  of  explicit  cloud  schemes. 


CHAPTER  I 


INTRODUCTION 


Clouds  are  among  the  most  important  regulators  of  the  weather  and  climate  in 
die  earth's  atmosphere.  Their  importance  for  the  atmospheric  circulation  has  been 
recognized  for  many  years.  However,  clouds  interact  with  atmospheric  processes  (i.e.. 
turbulence,  large  scale  circulation  and  radiation)  in  many  ways,  and  hence  their  role  is 
still  poorly  understood.  The  lack  of  knowledge  of  cloud  process  has  also  impaired 
cloud  parameterization  in  climate  and  forecast  models. 

During  the  last  decade  there  has  been  a  considerable  improvement  in  model 
representation  of  clouds.  Slingo  (1987)  improved  model-generated  cloud  fields  by 
linking  clouds  to  atmospheric  properties  such  as  relative  humidity,  vertical  velocity, 
and  static  stability.  This  scheme,  the  so-called  implicit  cloud  scheme,  is  successful  in 
producing  some  of  the  gross  features  of  global  cloudiness  but  it  lacks  a  sound  physical 
basis  and  does  not  represent  the  interaction  between  clouds  and  the  hydrological  cycle 
of  the  model.  However,  because  the  implicit  scheme  is  simple,  it  is  widely  used  in 
large-scale  models.  More  recently,  more  physically  realistic  schemes  (explicit  cloud 
schemes)  have  been  developed  (Sundqvist,  1978:  Lc  Treut  and  Li,  1989:  Smith,  1990; 
Tiedtke,  1991,  1993;  Mannoji,  1993).  With  these  schemes  cloud  liquid  water  content 
and  cloud  cover  are  predicted  through  a  prognostic  equation.  The  advantage  of  this 


1 


2 


•pimMch  is  that  it  represents  the  interaction  between  clouds  and  the  hydrological  cycle 
of  the  model. 

It  has  been  realized  (Tiedtke.  1993)  that  in  order  to  improve  the  model 
representation  of  clouds,  the  informadon  content  of  the  initial  state  must  also  be 
enhanced.  At  present,  operational  data  assimilation  systems  do  not  provide  a  robust 
initial  state  over  tropical  regions  due  to  the  sparsity  of  conventional  observations.  This 
causes  low  equilibration  of  the  humidity  field,  which  is  the  so-called  model  spin-up.  It 
reduces  forecast  skill  of  numerical  weather  predictions  (NWP)  at  the  lower  latitudes. 
Krishnamuiti  et  al.  (1984)  proposed  the  physical  initialization  scheme  which  improves 
the  initial  state  by  assimilating  observed  rain  rates.  Physical  initialization  restructures 
die  vertical  distribution  of  humidity  and  constrains  the  evolution  of  divergence 
distribution,  surface  pressure  fields  and  heating  from  the  condensation  process  to  be 
consistent  with  observed  rain  rates  and  model  physics.  In  impact  studies  of  physically 
initialized  data  sets,  Krishnamuiti  et  al.  (1991)  have  shown  an  overall  improvement  in 
die  distribution  of  the  rainfall  and  a  very  high  correlation  between  the  observed  and 
model  forecasted  rain  rates.  Mathur  et  al.  (1992)  also  have  shown  a  strong  positive 
impact  of  physical  initialization  on  tropical  rain  forecasts.  Krishnamuiti  et  al.  (1994) 
have  shown  that  the  correlation  for  1  day  forecast  with  the  physically  initialized  data  is 
0.55  which  exceeds  both  the  current  operational  skills  of  0.2  to  0.5  and  the  0.3 
correlation  of  persistence. 

It  is  our  contention  that  the  physically  initialized  data  set  with  its  enhanced 
information  will  have  strong  positive  impact  on  cloud  forecasts.  With  this  intention, 
we  designed  a  set  of  experiments  to  test  the  impact  of  physically  initialized  data  set  on 
two  cloud  parameterization  schemes,  namely,  the  implicit  cloud  scheme  and  explicit 
cloud  scheme.  Furthermore,  we  also  intended  to  illustrate  the  extreme  sensitivity  of 
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cloud  forecasts  to  the  initial  cloud  distribution.  Much  of  the  previous  related  work  has 
focused  on  improvements  to  the  parameterization  schemes.  However,  this  study 
emphasizes  the  fact  that  the  initial  data  is  equally  important  in  cloud  forecasts.  In  this 
regard,  efforts  at  European  Center  for  Medium-Range  Forecasts  (ECMWF)  are  already 
underway  to  incorporate  the  explicit  cloud  scheme  into  the  data  assimilation  system  so 
that  cloud  variables  are  carried  forward  in  time  (Tiedtke,  1993). 

The  thesis  has  been  divided  to  five  chapters.  The  model  is  described  in  Chapter 
IL  I^ysical  initialization  is  discussed  in  great  detail  in  Chapter  m.  The  numerical 
experiments  conducted  and  results  are  presented  in  Chapter  IV.  Hnally,  the 
conclusions  are  presented  in  Chapter  V. 


CHAPTER  II 


THE  FSU  GLOBAL  SPECTRAL  MODEL 


2.1  Overview 

Two  versions  of  the  Florida  State  University  Global  Spectral  Model  have  been 
used  in  this  study  on  cloud  prediction.  One  involves  an  implicit  cloud  scheme,  while 
the  other  utilizes  an  explicit  cloud  scheme.  Full  documentation  of  the  FSUGSM  is 
given  in  Krishnamurti  et  al.  (1991).  The  FSUGSM  includes  the  following  features: 

(a)  Independent  variables:  (x,  y,  or,  t). 

<b)  Dependent  variables:  vorticity,  divergence,  surface  pressure,  vertical 
velocity,  temperature  and  humidity. 

<c)  Horizontal  resolution:  Triangular  106  waves. 

(d)  Vertical  resolution:  14  layers  between  roughly  50  and  1000  mb. 

(e)  Semi-implicit  time  differencing  scheme. 

(f)  Envelope  orography  (Wallace  et  al,  1983). 

(g)  Centered  differences  in  the  vertical  for  ail  variables  except  humidity 
which  is  handled  by  an  upstream  differencing  scheme. 

(h)  Fourth-order  horizontal  diffusion  (Kanamitsu  et  al,  1983). 

(i)  Kuo-type  cumulus  parameterization  (Kuo  1965,  1974;  Krishnamurti  et  al, 
1983). 

(j)  Shallow  convection  (Tiedke,  1984). 

(k)  Dry  convective  adjustment. 
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(1)  Large  scale  condensation  (Kanamitsu,  1975). 

(nn)  Surface  fluxes  via  similarity  theory  (Businger  et  al,  1971). 

(n)  Vertical  distribution  of  fluxes  utilizing  diffusive  formulation  where  the 
exchange  coefficients  are  functions  of  the  Richardson  number  (Louis, 
1979). 

(o)  Long-  and  shortwave  radiative  fluxes  based  on  the  emissivity/absorptivity 
methods  (Chang.  1979). 

(p)  Diurnal  cycle  with  respect  to  the  radiative  processes. 

(q)  Implicit  cloud:  Parameterization  of  low,  middle  and  high  clouds  based  on 
thr^hold  relative  humidity  for  radiative  transfer  calculations 
(Kiishnamurti  et  aL,  1991). 

Explicit  cloud:  Parameterization  of  cloud  water  content  and  cloud  cover 
(Mannoji,  1993). 

(r)  Surface  energy  balance  coupled  to  the  similarity  theory  (Kiishnamurti  et 
al.,  1988). 

The  14  level,  7106  FSUGSM  used  in  the  study  has  a  spatial  resolution  of 
approximately  0.94  degrees  with  160  Gaussian  latitudes  and  384  Gaussian  longitudes. 
The  vertical  structure  is  defined  by  the  sigma  coordinate.  Figure  2.1  shows  a  schematic 
of  the  vertical  structure  of  the  14  level  model. 

Implicit  cloud  scheme  and  related  physical  processes  are  described  in  section  2.2. 
A  description  of  explicit  cloud  scheme  is  detailed  in  section  2.3.  Section  2.4  discusses 
the  radiative  transfer  processes  of  the  model,  which  is  a  common  feature  in  both 
schemes.  A  modification  to  parameterization  of  deep  cumulus  convection  in  the 
implicit  cloud  scheme  is  presented  in  section  2.5. 


2.2  Implicit  cloud  scheme 
2.2.1  Definition 

Clouds  are  classified  as  high,  middle,  and  low.  High  clouds  are  defined  to  occur 
between  400  and  100  mb,  middle  clouds  between  700  and  400  mb,  and  low  clouds 
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between  900  and  700  mb.  Clouds  are  not  defined  for  sigma  levels  below  low  cloud 
bases  or  above  high  cloud  tops. 

The  amount  of  a  cloud  type  in  an  individual  layer  is  assumed  to  be  a  function  of 
the  relative  humidity,  which  is  dependent  on  the  moisture,  temperature,  and  pressure. 

Clouds  are  assumed  to  be  present  when  the  mean  relative  humidity  (RH)  in  a  layer 
exceeds  the  threshold  value  (RHc).  The  cloud  amount  N  for  each  type  of  cloud  is 
defined  by  the  relation. 


N 


(h,  m  or  l) 


“  RH  -  RHc 
1  -  RHc 


RH>RHc 


(2.1) 


where  RHc  has  been  set  to  0.25,  0.40,  and  0.50  for  high,  medium,  and  low  clouds, 

respectively.  RH  is  the  pressure  weighted  mean  relative  humidity  in  a  layer.  The 

anaximum  possible  value  of  N  is  1.  The  cloud  amount  N  is  set  to  0,  when  RH  is  less 
than  RHc.  Total  cloud  amount  (N^  is  determined  from  a  random  overlap  of  the  cloud 
layers  in  a  column  as  given  by  N  =  1  -  (1-N,)(1-N  )(1-N  ).  Cloud  coverage  (CCn)  is 
partitioned  into  eight  sky  condition  categories  as  defined  in  Krishnamurti  et  al.  (1990): 


CCi  =  (l-NJ(l-Nj^)(l-N^) 

CC,  =  (  1  -  )  (  1  -  Njj ) 

CC,  =  N^N„(1.N„) 

CC,  =  N„(1-NJ(1.N„) 

CC,  =  N„N„(1.N^) 

CC,  =  N„(1.N^)(1-N„) 


clear  sky 

low  clouds  only 

low  and  middle  clouds 

low,  middle,  and  high  clouds 

middle  clouds  only 

middle  and  high  clouds 

high  clouds  only 


low  and  high  clouds 


CC.  =  N^N„(l-NJ 

This  sky  condition  is  used  in  the  computation  of  the  solar  radiation  absorbed  or 
scattered  by  clouds. 

2.2.2  Deep  cumulus  convection 

The  deep  cumulus  convection  of  the  FSUGSM  follows  Krishnamurti  et  al. 
(1993).  The  large  scale  supply  of  moisture  available  for  the  clouds  is  defined  by 


where, 

0)  =  large  scale  vertical  velocity 
q  =  large  scale  specific  humidity 
g  =  gravity 

Pt  =  pressure  at  cloud  top 
Pb  =  pressure  at  cloud  base 

In  this  study  an  additional  source  of  moisture  is  introduced  to  represent  the 
mesoscale  (subgrid  -  scale)  moisture  supply,  1^7).  The  total  moisture  supply  I  is  then; 

I  =  (1  +  Tj)  (2.3) 

Following  Kuo  (1974),  the  rainfall  rate  R  and  moistening  rate  M  may  be  expressed  as 
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M  =  Ib  =  I,  (1  +  Tj)  b  (2.5) 

where  b  is  a  moistening  parameter. 

FoUowing  Krishnamurti  et  al.  (1983),  the  supply  of  moisture  required  to  produce  cloud 
over  a  unit  area  is  as 


where 


Q  -  Qq  +  Qg 


(2.6) 


(2.7) 


(2.8) 


where  L  is  latent  heat  of  condensation,  and  Cp  is  specific  heat  of  dry  air  at  constant 
pressure.  The  term  is  the  moisture  required  tc  change  the  environmental  humidity 
distribution  from  q  to  q,  in  a  cloud  time  scale  (At)  of  30  minutes. 

Two  additional  proportionality  factors  are  now  defined: 


and 


1(1  -  b)  I  (1  +  tj)  (1  -  b) 

a  = - =  ± - 

Ib  ^ 


(2.9) 


(2.10) 


They  may  be  related  to  the  convective  heating  and  moistening  rate  as: 


—  +  V  •  V0  =  a. 

9t  • 


6.-6  de 

S - 2+  (0^ 

At  dp 


(2.11) 


^  +  V  •  Vq=:aq^-~-3^  (2.12) 

dt  At 


All  of  the  variables  to  determine  the  convective  heating  and  moistening  rates  are 
known  from  the  large-scale,  except  7}  and  b  (a  and  aq  are  determined  directly  from 
equations  (2.9)  and  (2.10)  once  7]  and  b  are  known).  Krishnamurti  et  al.  (1983) 
determined  7)  and  b  based  on  a  screening  multi-regression  aiudysis  of  GATE 
observations.  They  regressed  M/I  and  R/I  against  various  large-scale  variables  such  as 
vertical  velocity,  divergence,  vertical  shear  of  the  horiz  ontal  wind,  and  vorticity.  They 
found  that  the  vertically  averaged  large-scale  vertical  velocity  and  the  700  mb  relative 
vorticity  provided  the  best  estimates.  The  regression  equations  used  were 


=  ^iCtoo  +  bi©  +  Cl 


(2.13) 


=  ^2^700  +  ^20)  +  C2 


(2.14) 


where  ai  =  1.58  x  10^  sec,  bi  =  3.04  x  102  sec/mb,  C]  =  0.476,  a2  =  1.07  x  10^  sec,  b2 
=  1.07  X  102  sec/mb,  and  C2  =  0.87.  These  regression  constants  are  used  with  the  time 
dependent  values  of  the  vertically  averaged  vertical  velocity  and  the  700  mb  relative 
vorticity  to  produce  time  invarying  values  of  the  mesoscale  moisture  convergence  and 
moistening  parameters.  They  are  computed  from  the  following  equations: 
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(2.15) 

(2.16) 


2.2.3  Shallow  convection 

Shallow  convection  in  the  FSUGSM  follows  the  scheme  of  Tiedtke  (1984)  which 
assumes  a  simple  diffusive  fonn  for  the  parameterization  of  shallow  convection. 
Shallow  convection  is  invoked  for  a  conditionally  unstable,  undisturbed  situation  via 
K'theory.  This  parameterization  scheme  is  used  at  points  where  the  cloud  base  is 
below  700mb  and  the  air  is  conditionally  unstable.  The  cloud  top  may  not  extend 
above  700mb. 

The  moisture  and  the  temperature  equations  for  this  parameterization  are 


and 


5  =  5 

dx  Bp 


Jk 

di  Bp  * 


Bp 


(2.17) 

(2.18) 


where  Kq  =  =  Kg2p2,  K  =  25  luVs,  g  is  the  gravity,  and  p  is  the  density  of  air. 

Krishnamurti  et  al.  (1990)  showed  that  a  constant  value  of  K  does  not  produce  a 
smooth  continuity  in  the  fluxes  at  the  lifting  condensation  level  (LCL).  Instead,  the 
following  form  of  K  was  suggested 


(2.19) 
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where  Fq  (LCL)  is  the  planetary  boundary  layer  moisture  flux  at  the  LCL. 

A  similar  definition  of  was  given  for  the  same  reason.  This  parameterization 
of  shallow  convection  conserves  moisture.  The  total  convergence  of  moisture  flux  in 
the  cloud  layer  equals  Fq  (LCL).  These  shallow  clouds  are  nonprecipitating  and  do  not 
produce  saturation  on  the  large  scale. 


2.3  Explicit  cloud  scheme 

Explicit  cloud  scheme  adopted  in  FSUGSM  follows  Mannoji  (1993). 

The  basic  equations  for  the  local  change  in  the  cloud  water  mixing  ratio  (m)  and 
the  cloud  cover  (f)  prediction,  respectively,  are 


An  An^  An^. 

_=_s+_«  + 

dt  di  dl 


,  ^EC 

at  at 


(2.20) 


* 

_ _ _ ST  ^  CV  ^ _ W  _ rc 

A  A  at  A  at 


(2.21) 


where  subscripts  ST,  CV,  PP,  and  EC  denote  stratiform  cloud,  convective  cloud, 
precipitation,  and  ev{q)oration  of  cloud,  respectively. 

The  temperature  and  moisture  equations  can  then  be  expressed  by 


dT  A  rr^  L  f  ^ST  .  ^EC  ,  ' 

—  =  A(T)  -  —  - + - + - 

d  t  Cd  A  at  at 


^  =  A(q)  + 
dt 


'  ^ST  ^  ^EC  ^  ^ER  ' 

A  at  at 

^  j 


(2.22) 


(2.23) 
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where  A(T)  and  A(q)  represent  change  in  temperature  and  moisture  respectively.  These 
changes  are  produced  by  processes  that  are  not  associated  with  the  explicit  cloud 
scheme  such  as  advection.  and  turbulence,  etc.  Subscript  ER  denotes  evaporation  of 

^ST  *”sT  ^EC  ^EC 

rain.  In  practice  note  that - - - and - - - 

Bt  dt  dt  dt 

The  overall  processes  of  the  explicit  cloud  scheme  are  schematically  illustrated  in 
Rgrirc  2.2. 

2.3.1  Stratiform  clouds 

Following  Sundqvist  (1978),  the  formation  and  dissipation  of  stratiform  clouds 
are  calculated  from  changes  in  temperature  and  moisture  due  to  non-convecdve 
processes.  Unlike  Suhdqvist's  scheme  in  which  cloud  cover  is  diagnostically  calculated 
from  relative  humidity,  this  scheme  predicts  cloud  cover. 

in  the  computation  of  the  rate  of  change  of  cloud  water  and  cloud  cover  it  is 
assumed  that 

1)  over  the  grid  box  the  temperature  distribution  is  uniform  and  thus  the  saturated 
humidity  field  is  also  uniform, 

2)  in  the  overcast  area  q  =  qs  and  in  the  open  area  q  is  evenly  distributed  between  a 
minimum  value  and  q^, 

3)  cloud  water  is  uniformly  distributed  between  zero  and  a  maximum  value  over  the 
overcast  area, 

4)  changes  of  T  and  q  are  uniform,  and 

5)  cloud  forms  when  the  relative  humidity  exceeds  the  critical  values.  The  above 
mentioned  assumptions  are  applied  in  a  single  grid  box.  We  note  that  the  minimum 
value  of  moisture  and  the  maximum  value  of  cloud  water  should  be  derived  from 


known  model  variables. 


EVAPORATION 


Hguie  22. 


Overall  process  adopted  in  the  FSUGSM  explicit  cloud  scheme. 
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To  close  the  parameterization  of  the  cloud  water  and  cloud  cover  a  relationship 
between  the  cloud  variables  and  the  large-scale  variables  needs  to  be  prescribed. 
According  to  Mannoji  (1993)  this  relationship  is 

=-a*5p+ltai^  (2.24) 

RCpT2  p  RT^ 

or 

&j  -  Aq,  =  M  -  oAm  (2.25) 

where 

M  =  &j  +  Si*  ^  5r  (2.26) 

p  RT2 

and 

y=:  fLfaa  (2.27) 

RCpT2 

In  equation  (2.24)  Aq,  and  Am  denote,  respectively,  change  in  the  cloud  saturated 
humidi^  and  the  cloud  water  while  &},  dp,  and  ST  denote  changes  in  large-scale 
variables,  q,  p,  and  T,  respectively.  These  large-scale  changes  are  known.  The 
saturated  humidity,  q,,  is  also  a  known  large-scale  model  variable.  R  is  gas  constant 
and  e  is  equal  to  0.622. 

Stratiform  cloud  processes  consist  of  two  processes:  formation  and  dissipation. 
Equation  (2.24)  is  used  to  find  solutions  for  each  process.  Defining  Af  to  be  change  in 
the  cloud  cover,  Mannoji  (1993)  proposed  a  closure  for  Am  and  Af  based  on  a 
geometric  method  as  illustrated  in  Figures  2.3  and  2.4  for  the  formation  and  dissipation 
processes,  respectively. 


Formation  process 
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Figure  2.3  illustrates  that  when  water  vapor  increases  or  temperature  decreases 
(that  is,  M  is  greater  than  0)  water  vapor  condenres,  and  thus  it  turns  into  cloud  water  if 
saturated,  while  q,  may  increase  by  due  to  temperature  change.  This  means  that 
condensation  occurs  when  specific  humidity  exceeds  q,  lAq,. 

From  Figure  2.3  the  equations  for  Am  and  Af  are  readily  obtained  as 


Am  =  (5q  -  AqJ(f  +  |  Af)  (2.28) 

and 

Af  =  (5q  -  AqJ  ~  ^  (2.29) 

q*  *  Qmin 

where  Q„in  can  be  represented  by  the  known  variables  as 


or 


q,f  +  9iiULjLai(i  .f)=:q 
2 

1  -  f 


(2.30) 


(2.31) 


This  formulation  ensures  that  q^n  is  consistent  with  both  the  moisture  in  the  grid  box 
and  its  saturated  value. 

Using  equations  (2.25),  (2.28),  and  (2.29),  the  solutions  for  Af  and  Am  are 


and 


|^)<f+MA)  +  Mr  +  1  -  y  {)<f+2MA)+l}2  - 

^_(M  -  yAM)(l  -  f)2 
2(q  ,  -  q) 


472MA(f+MA) 


(2.32) 

(2.33) 


where 


C|s+'rAc]s''^ 

P  mox  "Qs 


q 


min - 


T-Sq 
— ^rAq 


s 


Hgure  2.3  Distribution  of  the  water  vapor  and  cloud  cover  within  a  grid  during 
the  stradfonn  cloud  formation.  Water  vi^r  distributions  before  (solid 
line)  and  after  (dariied  line)  condensation  occurs  are  shown.  Amount 
of  the  cloud  water  condensed  at  one  time-step  is  shaded. 


Figure  2.4  Distribution  of  the  water  vapor,  cloud  water  and  cloud  cover  in  a  grid 
during  the  stratiform  cloud  dissipation.  Shades  show  cloud  water 
amount  before  (heavily)  and  after  (lightly)  dissipation  occurs. 
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(2.34) 


4(q,  -  q) 

Given  die  large  scale  variations  Sp,  &],  and  ST,  Am  and  Sf  can  readily  be  calculated. 
Dissipation  ProccM 

Mannoji  (1993)  again  utilizes  the  geometric  method  to  achieve  the  closure  of  the 
dissipation  process  in  a  manner  similar  to  that  used  for  the  formation  process. 

Figure  2.4  shows  that  when  water  vapor  decreases  or  temperature  increases  (i.e. 
M  is  smaller  than  zero)  cloud  water  evaporates  and  thus  it  turns  into  water  vapor  if 
unsaturated,  while  q,  may  decrease  by  xAq^  due  to  temperature  change.  Evaporation, 
therefore,  occurs  when  specific  humidity  becomes  less  than  q^  +  xAq^  where  Aq,  is 
negative. 

The  equations  for  a  closure  of  Af  and  Am  are  then  given  by 


Af  _  f 
5q  -  Aq,  2m/f 
and 

2(m4-tAm)  /  (fftAf)  _  2m/f 
f  +  tAf  f 


(2.35) 


(2.36) 


where  2m/f  represents  a  maximum  value  of  cloud  water.  Equation  (2.36)  can  be 
written  as 


Ara  =  !S^(2  +  9M) 
f  f 


(2.37) 


Note  that  when  xAf/f  equals  -1,  the  cloud  disappears.  This  means  that  cloud  cover  and 
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cloud  water  becomes  zero  and  Am  becomes  equal  to  mAi/f  from  equation  (2.37). 
With  the  aid  of  equations  (2.2S).  (2.35)  and  (2.37).  the  solution  for  Af  is 


Af  = 


1_ 

Tf 


-(1  +  )f)  + 


Once  Af  is  determined,  Am  is  calculated  from  equation  (2.37). 


(2.38) 


2.3.2  Convective  clouds 

Total  condensation  in  a  column  and  time  tendency  of  T  and  q  at  each  level  are 
c<Mnputed  by  the  FSUGSM  cumulus  convection  scheme  described  earlier.  Two 
assumptions  are  required  to  calculate  the  cloud  water  source.  These  assumptions  are: 


1)  The  total  condensation  in  the  column,  P,  represents  the  vertically  integrated 
cloud  water  source  and  is  given  by 


P 


i  I 

dx  ^ 


(2.39) 


This  is  based  on  conservation  of  mass  of  water  in  the  column. 


2)  The  amount  of  cloud  water  at  each  level  is  proportional  to  a  vertical 
distribution  function,  D,  which  is  arbitrarily  chosen  to  have  a  linear  variation  from  1  at 
the  cloud  base  to  0.4  at  the  cloud  top.  Here,  cloud  top  and  cloud  base  are  provided  by 
the  implicit  cumulus  convection  scheme. 

For  a  complete  closure  of  the  convective  cloud  cover  formation,  the  cloud  cover 
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is  determined  following  Le  Trent  (1990).  He  assumes  that  the  uniformly  distributed 
total  water  within  the  grid  box  is  split  into  a  vapor  part  and  a  liquid  part,  provided  that 
the  total  water  and  its  maximum  and  minimum  values  are  given.  This  is  best  illustrated 
in  Figure  2.5. 

The  closure  equation  obtained  from  Rgure  2.5  is: 

=  9«  ~  9"»n  (2.40) 

f  1  -  f 

Combining  equations  (2.31)  and  (2.40),  the  equation  for  the  cloud  cover  due  to  the 
eonvecdve  cloud  formation  process  is  given  by 

1 

^  = -  '  (2.41) 

1  -f/  q,..:..,q 

^  m 


Assuming  that  the  cloud  cover  becomes  fcv  when  fcv  is  greater  than  the  cloud  cover  at 
the  previous  time  step,  f(t  -  1),  the  equation  can  be  written  by 


dt 


f(t  -  1) 


T 

0 


iffcv>f(t-  1) 

iffcy^f(t-  1) 


(2.42) 


2.3.3  Precipitation 

Mannoji  (1993)  used  the  Sundqvist  scheme  to  calculate  the  conversion  of  cloud 
water  to  raindrops,  cloud  ice  sedimentation,  evaporation  of  raindrops,  and  the  rate  of 
change  of  cloud  water  and  precipitation  fluxes. 

Since  the  model  does  not  separate  cloud  ice  and  cloud  liquid  water,  the 


18 


proportion  of  cloud  liquid  water  needs  to  be  prescribed.  Following  Sundqvist  (1989), 
this  is  diagnostically  prescribed  as: 


0 


I  1 


forT  <  Ti 

forTi  ^  T<T2  (2.43) 

for  T  ^  To 


where  Tj  and  T2  are  taken  to  be  233  K  and  273  K,  respectively.  These  are  used  for  a 
smooth  transition  between  frozen  and  liquid  states.  Cloud  in  the  precipitation  process 
is  divided  into  two  parts.  One  part  is  liquid  cloud  in  which  liquid  water  content  is 
given  by  Rm  and  cloud  cover  by  Rf.  The  other  part  involves  cloud  ice  in  which  ice 
content  is  represented  by  (1  -  R)m  and  cloud  cover  by  (1  -  R)f.  Therefore,  the  ratio  of 
cloud  water  to  ice  content  does  not  change  in  the  overcast  area. 

Cloud  to  raindrops 

Following  Sundqvist  (1989),  the  conversion  rate  of  cloud  water  to  rain,  Piiq,  is 
given  by 


Pu,  =  C,m{l-exp[-[!S^]  ]}  (2.44) 

where  C,  denotes  the  time  constant,  mg  a  critical  cloud  water.  In  practice,  this  equation 
is  modified  by  Sundqvist  (1989)  as 


where 


Puq  =  CtR  (m  -  fmcF)  F 


(2.45) 


Figure  2.5 


Distribution  of  the  water  vapor  and  cloud  water  to  determine  the  cloud 
cover  of  cumulus  cloud  in  a  grid.  Amount  of  the  cloud  water  is 
shaded. 
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>==>-“p[-[=x]  ] 

(2.46) 

Ct  =  Co  (1  +  Cl/pi)  {I+C2  7(268  -  T)  } 

(2.47) 

me  =  mo  /  [(1  +  CiTl^)  {1  +  C2  y  (268  -  T)  }] 

(2.48) 

where  Cq.  C],  C2,  and  mo  are  constants  listed  in  Table  2.1,  and  Pf  denotes  the  rate  of 
precipitation  which  will  be  shown  in  the  next  section. 


Table  2.1:  Parameters  used  in  parameterization  of  precipitation. 

(stratus/cumolus  cloud  for  Cq  tno) 


Co  lt>-3  / 10-*  s-i 

mo  5  X  1(H  /  3  X  1(H  kgAtg 


Cl  300 


(kgm-2s-i)’^^2 


C2  0.5 


K 


-1/2 


Cloud  ice  sedimentation 

The  rate  of  change  in  cloud  ice  due  to  flux  divergence  is  expressed  as 


(2.49) 
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where  the  flux  divergence  of  cloud  ice  is  given  by 


Fi  =  (1  -  R)  f(p,  mi  V)  (2.50) 

where  p,  is  the  density  of  air,  mj  is  cloud  ice  mixing  ratio,  and  v  is  the  fall  speed  of 
cloud  ice. 

Unlike  cloud  liquid  water,  the  fall  speed  of  cloud  ice  is  not  negligible. 
Heymsfield  (1977)  formulated  the  fall  speed  of  ice  from  observations  as 

v  =  3.326  (pmO*^’’  (2.51) 

The  conversion  rate  of  cloud  ice  to  precipitation  due  to  sedimentation  process  can 
tiien  be  expressed  as 


Pi  = 


S  for  S  >  0 
0  forS^O 


(2.52) 


Finally  the  rate  of  change  of  cloud  cover  due  to  precipitation  is  given  by 


(1-R)f/T  ifPiT>(l-R)ra 

£iX£=  . 

dt  0  ifPiT<(l-R)m 


(2.53) 


Evaporation  of  raindrops 

The  equation  for  evaporation  of  raindrops  is  expressed  in  Sundqvist  (1989)  as 
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at 


(2.54) 


Pr 


where  is  raindrop  evaporation  constant,  which  is  taken  to  be  10”^  (kg 

Rate  of  change  of  cloud  water  and  precipitation  flux 
The  rate  of  change  of  cloud  water  is  computed  from 


An, 


dt 


=  S-Pu,-Pi 


(2.55) 


and  the  flux  of  precipitation  is  given  by  the  equation 


-if 

'  «Jo 


Pliq  +  Pi 


ER 


at 


dp 


(2.56) 


Therefore,  total  precipitation  on  the  ground,  P,  is 


P  =  P,  +  Fi 


(2.57) 


2.3.4  Evaporation  of  cloud 

Evaporation  of  cloud  is  accelerated  by  the  turbulent  mixing  between  cloud  air 
and  unsaturated  environment  air.  This  is  taken  into  account  in  the  cloud  dissipation 
process  (Tiedtke,  1991)  as 


— ^  =  -K  f(q,-q) 

at 


(2.58) 
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^ere  cloud  evaporation  constant,  is  equal  to  l(Hs”'. 


(2.59) 


2.3.5  Radiation  for  cloud  liquid  and  ice 

The  predicted  cloud  water  content  is  employed  in  radiation  scheme  to  calculate 
the  optical  depth  and  emissivity.  Emissivity  is  computed  in  the  infrared  and  shortwave 
regions. 

In  the  infrared  region,  the  emission  function  is  computed  from  the  expression: 


s=  1  -  cxp(-kW)  (2.60) 

where  k  and  W  denote  the  absorption  coefficient  and  the  cloud  water  path,  respectively. 
The  absorption  coefficient  is  taken  to  be  0.14  for  cloud  liquid  and  0.05  for  cloud  ice 
(Stephens,  1978;  Starr  and  Cox,  1985). 

In  the  shortwave  region,  the  emission  function  is  calculated  from  the  equation 
given  by  Platt  and  Harshvardhan  (1988): 


1  -  exp(-0.757  Tj^)  (2.61) 

where  represents  the  optical  depth  in  the  shortwave  region.  Stephens  (1978) 
expressed  the  optical  depth  for  cloud  liquid  as 

^sw  “  ^  ^ 
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and  for  cloud  ice  as 

=  93.409  W  (2.63) 

sw 

2.3.6  Initial  condition  and  computational  procedure 

Initial  values  of  the  cloud  variables  are  based  on  the  procedure  given  by 

Pudyldewicz  et  d.  (1992).  Initially,  the  mean  relative  humidity  (RH)  is  taken  to  be  the 
relative  humidity  at  each  pressure  level  in  the  column.  Given  this  humidity  field,  the 
cloud  cover  is  diagnosed  from  equation  (2.1).  The  cloud  water  content  (CWC)  is 
ktitially  assigned  as 


CWC  =10  ”  inKg/m3  (2.64) 

Here,  the  maximum  value  of  cloud  water  content  is  limited  to  5  x  10~^  Kg/m^. 

To  avoid  unrealisdc  fall-outs  of  condensation,  discretion  in  computational 
procedure  is  advised.  The  sequence  of  processes  adopted  in  the  explicit  cloud 
parameterization  scheme  of  FSUGSM  is  as  follows:  1)  convective  cloud  formation,  2) 
stratiform  cloud  formation  and  dissipation,  3)  removal  of  supersaturation,  4) 
precipitation  process  and  raindrop  evaporation,  and  5)  evaporation  of  cloud  by  turbulent 
mixing. 

The  temperature  and  moisture  changes,  8l  and  5q,  are  computed  using  the 
large-scale  parameters,  dq  is  obtained  by  adding  all  changes  in  moisture  due  to  vertical 
advection,  turbulent  transfer,  and  adjustment  ST  is  calculated  in  a  similar  way,  but 
with  the  inclusion  of  changes  in  temperature  due  to  radiative  transfer. 

Adjustment  scheme  in  the  FSUGSM  explicit  cloud  scheme  is  described  in  Betts 
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(1973)  and  Smith  (1990).  During  the  changing  state  of  cloud  water,  the  liquid  water 
temperature,  Tu^,  and  the  total  water  content,  q,,  are  conserved.  These  are  defined  as 


and 


qt  =  q  +  m 


(2.65) 

(2.66) 


In  the  unstable  layers  the  lapse  rate  of  ,  is  set  equal  to  the  dry  adiabatic 

lapse  rate,  F,).  The  change  in  the  total  water  content  is  calculated  in  those  layers  in 
which  the  lapse  rate  has  been  changed  to  F^.  Since  the  total  water  content  is 
conserved,  th^  change  of  total  water  should  be  set  to  be  zero.  Temperature  and 
moisture  changed  by  the  adjustment  process  are  used  in  the  stratiform  cloud  process  to 
calculate  condensation  and/or  evaporation. 


2.4  Radiation 

The  radiation  parameterization  of  the  FSUGSM  consists  of  the  parameterization 
of  the  short-wave  and  long-wave  radiation.  These  are  briefly  described  in  the 
following  sections. 


2.4.1  Short-wave  radiation 

The  short-wave  radiation  parameterization  of  the  FSUGSM  is  based  on  the 
scheme  of  the  UCLA/GLAS  GCM  as  described  by  Davies  (1982)  and  Harshvardhan'er 
at.  (1987).  It  includes  a  parameterization  for  the  major  absorption  processes  in  the 
stratosphere,  the  troposphere,  and  at  the  earth’s  surface.  The  parameterization  utilizes 
the  water  vapor  distribution,  the  cloud  coverage,  the  zenith  angle  of  the  sun,  the  albedo 
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of  die  earth's  surface,  and  the  ozone  distribution. 

In  this  scheme  multiple  scattering  is  taken  into  account  Following  is  a  brief 
description  of  the  short-wave  radiation  parameterization.  The  heating  rate  at  pressure 
level  p  due  to  the  short-wave  radiative  flux  is  given  by 

S^(p)  -  SCp)-^)  (2.67) 

dt  Cp  dp 

I  + 

where  S  (p)  is  the  downward  solar  radiative  flux  at  pressure  level  p,  and  S  (p)  is  the 

upward  solar  radiative  flux  at  p  by  reflection  from  layers  below  p  and  the  ground. 

2.4. 1.1  Basic  equation  for  transfer  of  shortwave  radiation 

The  basic  equations  for  the  downward  and  upward  solar  flux  at  pressure  p  for  a 
nonscatteiing  atmosphere  are,  re^cdvely: 


S^(p)  =  /io  So  T(u(p)) 

(2.68) 

s^(p)  =  Po  So  Rg  T(u(p)*) 

(2.69) 

where 

Hq  =  cosine  of  the  zenith  angle  (6b) 

So  =  net  solar  radiation  at  the  top  of  the  atmosphere 
Rg  =  surface  albedo  integrated  over  the  entire  solar  spectrum 
^u(p})  =  mean  transmittance  integrated  over  the  entire  solar  spectrum 
u  =  optical  path  for  the  absorber 

u*  =  effective  optical  path  traversed  by  the  diffusively  reflected  radiation 


=  Muo  +  (uq  -  u)  M 


(2.70) 


where 


uo  s  total  optical  path  up  to  the  ground 
M  s  effective  magnification  factor  for  diffuse  radiation 


The  definitimi  of  t(u)  is  given  by 


where 


t(u) 


o 

=  —  f  exp  [-  M  f  k  du]  dv 

K  ’ 


(2.71) 


V  =  spectral  wavenumber 
Av  s  ^)ectral  width 

M  s  magnification  factor  accounting  for  the  slant  path  and  refraction  which 
is  given  by 


M  = 
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(1224  ^lo^  +  1) 


s  monochromatic  mass  absorption  coefficient 


(2.72) 


Different  atmospheric  constituents  (water  vapor,  O3,  etc.)  have  different  absorption 
properties.  Since  the  heating  rate  by  the  short-wave  radiation  can  be  expressed  in  terms 
of  absorption  functions,  it  is  important  to  determine  the  absorption  functions. 


2.4. 1.2  W ater  vapor  absorption 

Absorption  by  water  vapor  is  the  major  source  of  solar  radiative  heating  in  the 
atmosphere  (Lacis  and  Hansen,  1974).  For  clear  skies  the  absorption  due  to  water 
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vapor  is  modeled  based  on  the  k-distribution  approximation  (Stephens,  1984)  in  which 
the  absorptivity,  Awv(y),  is  expressed  as 

Awv(y)=l-|  p(k)e-kydk  (2.73) 


where, 

p(k)  =  probability  distribution  for  the  absorption  coefficient 
p(k)dk  =  fraction  of  the  incident  flux  that  is  associated  with  an  absorption 
coefficient  between  k  and  k  +  dk 
y  =  effective  optical  depth  of  water  vapor 

The  effective  optical  depth  of  water  vapor  is  given  by 


where, 

M  =  magnification  factor 

g  =  gravity 

q  =  specific  humidity 

Po  =  reference  pressure  =  1013.25  mb 

To  =  reference  temperature  =  273.16  K 


(2.74) 


For  the  upward  beam  in  a  clear  sky,  the  effective  optical  depth  is  given  by 
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The  net  absorption  by  water  vapor  in  a  clear  sky  in  the  i***  layer  is  then: 

Aiwv  ~  AkjSo  ■|Ay^(yj^j)  -  Ay^(yi)  +  R^Ot)  [A^(yj^j)]|  (2.76) 

where 

So  =  solar  flux  at  the  top  of  the  atmosphere 

Rf  (/i)  =  ground  albedo 

For  cloudy  skies  multiple  scattering  is  the  principal  factor  affecting  the 
absorption  by  water  vapor.  The  parameterization  of  cloudy  sky  water  vapor  absorption 
is  also  based  on  the  discrete  probability  distribution  which  is  the  k-distribution  method. 
Tlie  ccmiplete  radiative  transfer  problem  including  multiple  scattering  is  solved  as 
discussed  by  Davies  (1982).  Apart  from  the  absorption,  the  reflection  and  transmission 
characteristics  of  each  atmospheric  layer  need  to  be  determined,  given  the  total  optical 
thickness  S,  the  single  scattering  albedo  a  ,  the  asymmetry  scattering  factor  g,  and  the 
solar  zenith  angle  6.  The  solution  is  done  through  a  two-stream  approximation.  From 
among  the  various  two-stream  approximation  available,  the  delta-Eddington 
approximation  given  by  Joseph  (1976)  is  used.  In  this  approximation,  the  transfer  of 
the  direct  solar  beam  radiation  through  a  single  homogeneous  layer  is  assumed  to  give 
rise  to  a  diffuse  radiation  at  the  top  of  the  layer  (the  reflected  part,  UO)  and  a  diffuse 
radiation  at  the  bottom  of  the  layer  (the  transmitted  part,  DO).  The  upward  diffuse 
radiation  UO  and  the  downward  diffuse  radiation  DO  are  given  as  a  function  of  S',  co' 
and  g'  where 
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5'  =  5(1  -  im 

S'  =  5)(1  - 12)  /  (1  -  pQ» 

i'  =i/(i  +  i) 


(2.77) 

(2.78) 

(2.79) 


All  these  parameters  were  transformed  based  on  the  asymmetry  factor  (g)  to  normalize 
the  delta*Eddington  phase  function.  The  reflection  R  and  transmission  T  functions  for 
each  homogeneous  layer  are  approximated  by  the  two-stream  approximation  of  Sagan 
and  Pollack  (1967); 


p  _(M-t-l)(u-l)(et  -  e-t) 
(u+l)2e‘  -  (u-l)2e-t 


(u+DV  -  (u-l)V* 

where 


(2.80) 

(2.81) 

(2.82) 


t=[3(l-S)(l-iS)]i5 


(2.83) 


Using  these  equations,  the  reflection  and  transmission  functions  and  the  upwelling  and 
downwelling  radiation  due  to  the  direct  solar  beam  can  be  calculated  at  all  levels. 

The  direct-beam  radiation  reflected  and  transmitted  by  a  layer  may  be  considered 
as  sources  of  diffuse  radiation  for  the  layers  above  and  below  it.  The  upwelling  and 
downwelling  radiation  at  each  layer,  as  given  by  Harshvardhan  et  al.  (1987),  is: 


UOi  =  ^SiRi()io) 


i  =  1, . . . ,  n+1 


(2.84) 
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DOi  =  ^iTiOio)  i  =  l . n  (2.85) 

where 

UOj  =  upwelling  radiation  through  level  i 
DO,  =  downwelling  radiation  through  level  i 
Ri  =  reflectivity  of  layer  i  to  unit  diffuse  radiation 
Ti  =  transmissivity  of  layer  i  to  unit  diffuse  radiation 
/ioS,  s  direct  solar  flux  incident  on  level  i 

The  downwelling  diffuse  radiation  at  level  i+1  due  to  all  diffuse  radiation  that  has  not 
creased  level  i  is 


Dli  =  (Dli  +  UO,  CRi)  MiTi  +  DOi  i  =  1. . .  .  ,  n+1  (2.86) 

where 

CRj  =  composite  reflection  coefficient  of  the  atmosphere  above  i 
M,  =  magniHcation  factor  for  multiple  reflections 

Here,  CR,  is  recursively  obtained  by 

CR,  =  Rj-i  +  Tj_i  CRj-i  Ti_i  Mi_i  i  =  2, ....  n+1  (2.87) 


Mj  is  given  by 


Mi  =  1  /  (1-CRiRi) 


i  =  1, .  . .  ,  n+1 


(2.88) 
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By  definition 


CRi=0  and  DOi=0  (2.89) 

The  upwelling  diffuse  radiation  at  i  due  to  all  diffuse  radiation  that  has  not  previously 
crossed  level  i+1  is 


Uli  =  (DljRi  +  UOi)  Mi  i=l,...,n+l  (2.90) 

Once  Dlj  and  Uli  are  found  for  all  i,  the  total  radiation  crossing  each  level  may  be 
recursively  calculated  by 

Ui  -  Uli  +  Ui+i  Ti  Mi  i  =  n . 1  (2.91) 

Di  =  Dli  +  Ui  CRi  +  /ioSi  i  =  n+1 . 1  (2.92) 

To  obtain  the  broadband  values,  each  of  the  spectral  values  calculated  above  are 
summed.  The  absorbed  radiation  within  each  layer  is  then  found  directly  from  the  net 
flux  convergence  of  solar  radiation. 

2.4. 1.3  Ozone  absorption 

According  to  Lacis  and  Hansen  (1974)  the  major  source  of  stratospheric  heating 
is  the  absorption  of  solar  radiation  by  ozone.  The  absorption  due  to  ozone  occurs  at 
wavelengths  where  Rayleigh  scattering  is  significant.  Most  of  the  Rayleigh  scattering 
takes  place  below  the  ozone  layer,  so  the  lower  atmosphere  acts  mainly  as  a  reflecting 
layer. 
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The  absorption  due  to  ozone  is  different  in  the  weak  visible  bands  (Chapius 
band)  versus  the  ultraviolet  bands  (Hartley  and  Huggins  bands).  The  percentage  of 
solar  flux  absorbed  as  function  of  ozone  amount  for  the  two  spectral  regions  are  given 
in  Lacis  and  Hansen  (1974).  They  showed  that  the  ultraviolet  bands  are  practically 
saturated  at  0.5  cm  of  ozone,  while  absorption  in  the  visible  bands  remain  nearly 
proportional  to  the  ozone  amount  The  total  fractional  absorption  is  expressed  as 


A(X)  =  Auv(X)  +  Avu(X) 


(2.93) 


where  X  is  effective  ozone  amount  in  cm  (STP),  Auv(X)  is  the  fraction  of  incident  solar 
flux  that  is  absorbed  in  the  ultraviolet  bands,  and  Avis(X)  is  the  fraction  of  incident 
solar  flux  that  is  absorbed  in  the  visible  bands. 

Lacis  and  Hansen  (1974)  fitted  the  frequency-integrated  absorption  curves  for 
ozone  to  analytical  expressions: 


A„v(X)  = 


1.082X  , 

(1+138. 6X)®*°* 


0.0658X 
1  +  (103. 6X)^ 


(2.94) 


Avi,(X)  = 


0.02118X _ 

1  +  0.042X  +  0 . 000323X2 


(2.95) 


The  ozone  amount  traversed  by  the  direct  solar  beam  in  reaching  the  i***  layer  of  the 
atmosphere  is 


where 


Xj  =  UiM 


(2.96) 


Uj  =  amount  of  ozone  in  a  vertical  column  above  the  i'^  layer 
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M  s  magnification  factor  accounting  for  the  slant  path  and  refraction 

The  ozone  path  traversed  by  the  diffuse  radiation  illuminating  the  i<^  layer  from  below 
is 


X*=u,M+ 1.9(u,-Ui)  (2.97) 

where  Ut  is  the  total  ozone  amount  in  a  vertical  path  above  the  main  reflecting  layer 
(which  is  the  ground  for  clear  skies  or  the  cloud  top  for  cloudy  skies).  The  total 
absorption  of  short-wave  radiation  in  the  i<>>  layer  due  to  ozone  is  then: 

Ai  =  /ioSo  {A(Xi+i)  -  A(Xi)}  +  ^ioSoR(A))  {a(X.)  -  A(X.^,)}  (2.98) 

where  R(/io)  is  the  albedo  of  the  reflecting  region  including  the  effective  albedo  of  the 
lower  atmosphere  and  the  ground  reflectivity. 

2.4. 1.4  Earth's  surface  absorption 

The  surface  absorption  is  divided  into  two  parts.  One  portion  concerns  the 
wavelength  regions  where  the  absorption  coefficient  of  water  vapor  is  significant  The 
other  part  covers  absorption  in  the  remaining  wavelengths.  Therefore,  the  total  surface 
absorption  is  the  sum  of  the  absorption  in  the  spectral  regions  associated  with 
significant  water  vapor  absorption  and  the  spectral  regions  of  negligible  water 
absorption.  That  is 

Ag  =  Agi  +  Ag2 


(2.99) 
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where 

Ag  s  total  surface  absorption 

Agi  =  surface  absorption  in  the  spectral  region  associated  with  significant  water 
vapor  absorption  (calculated  differently  for  clear  and  cloudy  skies) 

Ag2  =  surface  absorption  in  the  spectral  regions  of  negligible  water  vapor 
absorption 

For  clear  skies. 


A,i  =  /ioSo  [1  -  R,(A))]  f  P(k)e-ky  dk  (2.100) 

’'o 

=  AZoSo  [1  -  Ro(Po)]10.647  -  Afu^M)]  (2.101) 

where 

Rg  s  ground  albedo 

s=  albedo  of  the  ground  and  the  atmosphere  above 

For  cloudy  skies,  A^j  is  obtained  by  multiplying  the  transmission  function  at  the 

ground  by  [1  -  Rg(Aio)].  Ag2  is  computed  using  equation  (2.99),  where  is 

obtained  accounting  for  the  optical  depth  for  clouds,  5c,  and  the  asymmetry  scattering 
factor  for  clouds,  gc- 

2.4.2  Long-wave  radiation 

The  long-wave  radiation  scheme  of  the  model  is  also  based  or  the  UCLA/GLAS 
CK)M  scheme  and  is  described  by  Harshvardhan  and  Corsetti  (1984).  It  includes  a 
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panuneterization  of  the  transmission  functions  for  carbon  dioxide  and  water  vapor  using 
the  band- wing  scaling  approximation  developed  by  Chou  and  Peng  (1983)  and  Chou 
(1984).  It  also  includes  the  computation  of  the  transmission  functions  for  ozone  based 
on  the  formulation  given  by  Rodgers  (1968). 

The  radiative  cooling  rate  due  to  the  long-wave  radiative  flux  is  given  by 


dT  _  g  d  (F4>(p)  -  FT(p)) 
d  t  Cp  dp 


(2.102) 


where  F4'(p)  and  FT(p)  are  the  long- wave  downward  and  upward  fluxes  at  pressure 
level  p,  respectively.  This  is  the  final  output  of  the  radiation  scheme. 


2.4.2. 1  Basic  equation  for  transfer  of  longwave  radiation 

The  basic  equations  for  the  clear-sky  longwave  downward  (F^ir)  and  upward 
(fIit)  fluxes  integrated  over  the  spectral  range  Av  are 


(2.103) 


(2.104) 


where 

B^(T,)  =  blackbody  flux  at  the  surface  temperature  T* 
Ps  =  pressure  at  the  bottom  of  the  atmosphere 
Pt  =  pressure  at  the  top  of  the  atmosphere 
T(p')  =  air  temperature  at  pressure  p' 
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T  (p.p')  =  diffuse  transmittance  between  levels  p  and  p' 


Integrating  by  parts  above  equations  (2.103)  and  (2.104)  gives 


T(pt) 

fL(p)  =  B(T(p))  -  G(p.pc.T(Pt))  +  f  dT(p')  (2.105) 

T(p) 


T<P,) 


fL(P)  =  B(T(p))  +  G(p.p..T,)  -  G(p.p,.T(p.))  +  [  ^(P»P'»'^iP-l^T(pO  (2.106) 

j  dfi 

T(P) 


where 


B(T(p))=f  B^(T)dv 


(2.107) 


G(p,p',T)=f  T(p.pOB^(T)dv 


(2.108) 


r  ,(p,p-)^dv 

rTT  J  Am  ^ 


(2.109) 


T(p,)  =  air  temperature  at  surface  pressure  p* 


If  the  spectral  width  Av  of  the  band  is  sufficiently  narrow,  B^(T)  may  be  replaced  by  a 
mean  value  of  the  Planck  function.  The  equation  for  G  becomes 

G(p.pM)  =  f  T  (p,p')  dv  (2.110) 

Av  Jav 

In  order  to  evaluate  the  G  function  for  the  fluxes,  it  is  necessary  to  determine  the 
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diffuse  transmittance. 

Chou  and  Peng  (1983)  and  Chou  (1984)  defined  the  molecular  line  absorption 
and  the  e-type  absorption  (continuum  absorption)  for  the  computation  of  the 
diffiise-transmittance  function  due  to  water  vapor  and  CO2.  The  spectral  ranges  for 
both  the  molecular  line  absorption  and  the  e-type  absorption  ate  classified  as  the 
band-centers  (i.e.,  regions  of  strong  absorption)  and  the  band- wings  (i.e.,  regions  of 
moderate  absorption).  Following  Chou  and  Peng  (1983)  and  Chou  (1984),  broadband 
transmission  functions  for  each  of  the  absorbers  is  derived  by  averaging  the 
diffuse-transmittance  function  associated  with  molecular  line  absorption  and  the 
diffuse-transmittance  associated  with  e-type  absorption  over  wide  spectral  intervals. 
The  Planck-weighted  broadband  transmission  functions  are  used  for  the  water-vapor 
bands,  while  for  the  CO2  bands,  the  broadband  transmission  functions  are  computed  by 
integrating  the  diffuse-transmittance  over  the  entire  spectral  interval.  As  a  further 
amplification,  Chou  and  Peng  (1983)  and  Chou  (1984)  fitted  the  broadband 
transmission  function  for  water  vr^r  and  CO2  to  analytical  functions. 

The  broadband  transmission  functions  for  O3  is  based  on  the  formulation  given 
by  Rodgers  (1968).  Rodgers  fitted  the  integrated  absorption  of  9.6  /im  O3  band  to  an 
analytical  form  with  the  Lorentz  line-shape  transmission  function  for  the  random 
model.  With  these  simplifications,  the  regression  relations  for  the  broadband 
transmission  function  for  O3  are  obtained. 

2.4.2.2  Water '  .’por  bands 

According  to  Chou  and  Peng  (1983),  the  diffuse  transmittance  associated  with 
molecular  line  absorption  between  pressure  levels  p  and  p '  at  wavenumber  v  is 
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T^(p.p')  =  ^Jy*)  =  if  cxp[-k(ppT,)  w(p,p'y/ii  ^  d/i 
*^0 

where 

k  s  molecular  line  absoiption  coefficient 

ft  s  cosine  of  the  zenith  angle 

p,  =  reference  pressure  ftom  Chou  (1984) 

Tr  =  reference  temperature  from  Chou  (1984) 
w  =  scaled  water  vapor  amount  and  is  given  by 


P2 

*(PI.|)2)  =  J I  € )  expircr  -  T,))  q(p)  dp 


where 


r  s  factor  from  Chou  (1980) 
q  s  water  vapor  mixing  ratio. 


The  diffuse  transmittance  associated  with  e-type  absorption  is 


where 


T^CP.P')  =  T^(u)  =  2  «P[-<y^(To)  u(p,p')//i]  M  d/i 

”^0 


a  =  e-typc  absorption  coefficient 
u  =  scaled  water  vapor  amount  and  is  given  by 


«(Pi.P2)  =  If  e(p)  exp[1800  ( -  f  )]  q(P)  dp 

*  Jpj  0 

where 

e(p)  =  water  vapor  pressure  in  atmospheres 


(2.111) 


(2.112) 


(2.113) 


(2.114) 
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Therefore,  the  "Planck  weighted"  transmission  function  is 

x(w.u;T)  *  --i-  f  B  m  T  (w)  T  (u)  dv  (2.1 15) 

B(D  J 

Since  r^(w,u;T)  is  a  slowly  varying  function  of  temperature,  it  can  be  fitted  by  the 
quadratic  function  (Chou.  1984).  The  regression  equation  for  broad  transmission 
function  for  the  band-wings  is  then: 

t^(w,u;T)  =  t(w,u;250)  [1  +  a(w,u)  (T  -  250)  +  ^(w.u)  (T  -  250)^]  (2.1 16) 


where 

a(w,u)  s  regression  coefficient,  given  by  Chou  (1984) 

/3(w,u)  =  regression  coefficient,  given  by  Chou  (1984) 

«(w,u,250)  s  standarxlized  transmission  function,  given  by  Chou  (1984) 

The  regression  equation  for  t^(w;T)  for  the  band-centers  is  the  same  as  for  the 
band-wings  but  for  t(w;250),  a(w)  and  ^w). 

The  G  function  for  the  water- vapor  band  can  be  evaluated  as 

where 

=  spectrally  integrated  Planck  blackbody  functions  over  the  band-wings 

WVw 

B  =  spectrally  integrated  Planck  blackbody  functions  over  the  band-centers 

WVe 


As  mentioned  earlier,  the  band-center  absorption  due  to  water  vapor  can  be  neglected 
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it  is  dominated  by  molecular  line  absorption  and  the  effect  of  e-type  absorption 

is  small. 

For  the  IS  region,  Chou  (1984)  fitted  the  diffuse  transmission  function 
integrated  over  the  entire  band.  The  diffuse  transmission  functions  for  water  vapor 
molecular  line  and  e-type  absorption  are,  respectively, 

T(w)  =  exp  [-6.7w/(l  +  Ibw***)]  (2. 1 1 8) 

t(u)  =  exp  [-27u®  ”]  (2.119) 

The  G  function  is  given  by 

=  (2.120) 

For  the  9.6  /zm  region,  the  molecular  line  absorption  is  weak  and  thus  only  the 
e-type  absorption  is  considered.  Chou  (1984)  also  fitted  the  diffuse  transmission 
associated  with  the  e-type  absorption  in  this  region  as 

t(u)  =  exp  [-9.79u]  (2.121) 


The  G  function  is  given  by 

°WV,.6  = 


(2.122) 


The  G  functions  for  water  vapor  in  both  the  spectral  regions  represent  part  of  the  total 
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G  function  required  to  calculate  the  upward  and  downward  longwave  radiation  fluxes. 

2.4.2.3  CO2  bands 

The  diffuse  transmittance  for  CO2  is  obtained  by  using  equation  (2.1 1 1).  That  is 


=  T^(wc)  =  2[  exp  (-k^(p^T,)  Wc(p.p'yMl  d/i  (2.123) 


where 


WcCp.p')  =  scaled  CO2  amount  and  is  ^ven  by 


P2 

Wc(Pi.P2)  =  i[  c(p)f(p.T(p))dp  ,  (2.124) 

where 

c(p)  a  CO2  concentration 

Equation  (2.123)  treats  an  inhomogeneous  atmosphere  as  a  homogeneous  atmosphere 
with  a  constant  pressure  and  temperature  of  Pr  and  .  This  is  possible  because 
equation  (2.124)  scales  the  CO2  concentration  to  simulate  the  absorption  in  an 
inhomogeneous  atmosphere. 

Chou  and  Peng  (1983)  fitted  the  spectrally  averaged  diffuse  transmittance  to 
analytical  function.  Their  fitting  function  was 

z(Wc)  =  exp  [-awc/(l  +  bw"  )]  (1  ’ 

where  a,  b,  and  n  are  constants  chosen  for  individual  spectral  intervals.  Parameter  a  is 
taken  to  be  3.1  for  the  band-centers,  0.08  for  the  narrow  band-wings,  and  0.04  for  the 
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wide  band-wings.  The  constant  b  is  taken  to  be  IS.l  for  the  band-centers,  and  0.9  for 
the  band-wings.  Lastly,  the  parameter  n  is  taken  to  be  0.56  for  the  band-centers,  0.55 
for  the  narrow  band-wings,  and  0.57  for  the  wide  band-wings. 

In  order  to  account  for  radiative  cooling  due  to  CO2  in  the  stratosphere,  Chou 
and  Peng  (1983)  used  the  scaling  function  given  by 

f 

f(P.T)=[£J  R(T,T,)  (2.126) 

where. 

m  =  parameter  for  correcting  the  error  arising  from  the  assumption  of  linear 
dependence  of  the  absorption  coefficient  on  pressure. 

RCT.Tj)  =  temperature  scaling  factor  and  is  given  by 

R(T,T,)  =  exp(r(T-Tr)]  -  (2.127) 

where 

r  =  0.0089  for  the  band-centers  and  0.025  for  the  band-wings 
Tr  =  240  K  for  both  the  band-centers  and  the  band-wings 

Since  the  pressure  scaling  in  equation  (2.126)  assumes  that  the  absorption 
coefficient  follows  the  Lorentz  line-shape  transmission  function  (equation  (2.129)),  this 
assumption  is  not  valid  at  low  pressures  where  broadening  of  absorption  lines  due  to 
Doppler  shift  is  important.  In  the  15  fjm  spectral  region  the  height  where  line 
broadening  due  to  molecular  collision  and  Doppler  shift  are  equally  important  is 
approximately  10  mb.  For  a  Doppler  line  function  the  absorption  is  independent  of 
pressure.  To  account  for  the  Doppler  effect,  we  define  a  critical  pressure  level,  pg.  At 
pressures  higher  than  this  critical  pressure,  the  absorption  coefficient  is  independent  of 


43 


imssuie.  We  may  now  rewrite  equation  (2.126)  as 

fCp.T)  =  [  ^  ]  RCr.T,)  .  for  p  >  Pc  (2. 128) 

vriiere  m  is  0.85  for  the  band-centers  and  the  narrow  band-wings,  O.SO  for  the  wide 
band-wings,  and  Pc  is  1  mb  for  both  the  band-centers  and  band-wings  (Chou  and  Peng, 
1983). 

2.4.2.4  Ozone  bands 

Rodgers  (1968)  defines  a  transmission  function  for  a  Lorentz  line  shape  as 

=  (2.129) 

where 

a  =  line  width  at  one  atmosphere  in  cm~^ 

S  =  mean  spectral  interval  in  cm~^ 

kj  =  line  strength  in  cm/g 

Wo  =  ozone  concentration  in  g/cm^ 

With  further  simplification.  Harsh vardhan  and  Corsctti  (1984)  fitted  the 
transmission  function  to  an  analytical  function  as 


i 

T(wo)  =  1  -  0.677  [1  -  exp  {-4.398p(l  +  79.4  (-0)  )  -  1}] 

P 


(2.130) 


where 
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P2 


p  =  effective  broadening  pressure  is  given  by 


P  = 


.P2 


1 

O3  pdp 
Pi 

1013.25 

1 

P2 

O3  dp 
Pi 

(2.131) 


(2.132) 


where  O3  is  the  ozone  mixing  ratio. 

In  summary,  broadband  transmission  functions  in  this  model  are  separately 
calculated  for  three  bands:  water  vapor,  15-;im,  and  9.6-;im  bands.  Water  vapor  band 
is  partitioned  into  two  bands.  One  partition  is  for  the  band-centers  in  which  the 
molecular  absorption  is  much  stronger  than  e-type  absorption,  so  that  e-type  absorption 
can  be  neglected  for  the  band-centers.  The  other  is  for  the  band-wings  where  both 
types  of  absorptions  are  included.  In  the  15-;im  band,  CO2  absorption,  water  vapor 
molecular  absorption,  and  water  vapor  e-type  absorption  are  present  In  the  9.6-/im 
band,  water  vapor  molecular  absorption  is  neglected.  Thus  the  total  absorption  is  due 
to  ozone  and  water  vapor  e-type  absorption.  Given  the  broadband  transmission 
functions,  the  spectrally  integrated  Planck  blackbody  function  over  each  band  can  be 
obtained.  Finally,  the  spectrally  integrated  Planck  blackbody  function  for  the  upward 
and  downward  long-wave  radiation  fluxes  is  obtained  by  adding  functions  for  each 
band.  Formally,  this  is  given  by 


B(T(p))  =  mp))  +  B^^(T(p))  +  B,^^a(p))  +  B,^(T(p))  (2.133) 
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Similarly,  the  total  G  function  is  given  as  a  sum  of  G  function  for  each  band.  That  is 


G(p,p'.T)  =  G  +  G,,_  +  G, 


lS>im 


9.6^ 


(2.134) 


A  similar  expression  for  dG/dT  can  also  be  obtained. 

To  account  for  the  effect  of  cloudy  skies,  the  probability  for  a  clear  line  of  sight 
is  defined.  The  equations  for  cloudy  skies  are  given  by 


T(p,) 

fL(P)  =  B(T(p))  -  C(p.p,)  G(p.pt.T(p,))  +  J  C(p.p') 

T(p) 


^G(p.p  .T(p  )  j j 


(2.135) 


T(p,) 

fL(p)  =  B(T(p))-C(p.p^  [G(p.p,.T^-G(p.Ps,T(p,))]  +  f  C  p')  ^G<2£liI(E:)<lT(p') 

^(p)  ^ 

a  136) 


where  C(p,p')  is  the  probability  for  clear  line  of  sight  from  p  to  p'. 

This  probability  depends  on  the  kind  of  cloud  overlap  present  between  two 
pressure  levels.  In  this  model  a  random  overlap  is  used  to  compute  the  probability  for 
clear  line  of  sight  It  is  equal  to  the  product  of  all  fractional  cloud  amounts  for  all 
levels  that  have  clouds  (i.e.,  N(3)  x  N(4)  x....).  It  should  be  noted  that  the  minimum 
clear  line  of  sight  for  a  random  overlap  atmosphere  is  (1  -  N(3))  x  (1  -  N(4))  x.... 

2.5  Modification  to  cumulus  convection 

An  examination  of  model  forecasted  cloud  amount  and  type  showed  dissipation 
of  clouds  with  the  time  (Krishnamurti  et  al.,  1988).  The  reason  for  this  was  found  to 
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be  the  drying  of  atmosphere  in  the  upper  layers.  A  careful  analysis  of  model  output 
showed  that  the  drying  was  occurring  primarily  in  regions  of  tropical  convection. 

The  FSUGSM  has  a  modified  Kuo  cumulus  parameterization  scheme  for  deep 
convection  (detailed  in  section  2.2.2).  The  drying  in  the  upper  layers  of  the  model 
atmosphere  was  due  to  excessive  removal  of  the  moisture  needed  for  precipitation  and 
moistening  processes  by  the  convective  scheme.  The  scheme  was  modified  to  remove 
more  moisture  from  lower  levels  with  higher  specific  humidity  than  from  the  upper 
levels.  It  was  found  that  by  scaling  the  moistening  term  in  each  layer  by  a  function 
proportional  to  the  undesirable  dissipation  of  clouds  was  removed.  Thus  the 
change  in  specific  humidity  at  any  level  due  to  convective  process  is  given  by 


(Aq),  =  a,i3^) 

Qc 


(2.137) 


where  k  represents  a  model  level,  and  is  a  column  integral  of  specific  humidity 
given  by 


''Pt 


dp 


(2.138) 


The  remaining  symbols  are  explained  in  section  2.2.2.  Therefore,  the  total  moisture 
change  removed,  Qqr,  is  given  by 


(2.139) 


or 


This  fonnulation  results  in  the  moisture  needed  for  moistening  being  removed  from  the 
taiyers  with  higher  water  content. 

The  above  modification  to  the  cumulus  convection  parameterization  was  tested 
using  the  T42  FSUGSM.  The  modified  scheme  achieves  the  desired  end  in  that  the 
resulting  globally  averaged  total  cloud  amount  remains  steady.  This  is  in  accord  with 
what  would  be  expected  in  the  real  atmosphere. 


CHAPTER  m 


PHYSICAL  INITIALIZATION 


3.1  Introduction 

Although  the  current  data  coverage  in  the  operational  World  Weather  Watch 
(WWW)  includes  surface  based  networks  of  radiosonde,  rawinsonde  and  pilot  balloon 
distributions,  diverse  asynoptic  surface  and  space  based  system  such  as  commercial 
•hips  and  aircrafts  observations,  soundings  from  the  polar  orbiting  satellites,  and 
cloud-tracked  winds  from  the  geostationary  satellites,  there  still  exist  large  data  gaps 
over  both  continental  and  oceanic  areas  of  the  tropics  as  shown  in  ICrishnamurti  et  al. 
(1991).  The  initial  conditions  used  in  NWP  over  these  areas  are  largely  derived  from 
the  first  guess  provided  by  a  short  forecast  from  a  previous  cycle  of  the  data 
as^ilation  system.  However,  due  to  model  imperfections  the  quality  of  the  analyses 
over  data  gap  regions  are,  therefore,  likely  to  be  lower  than  over  data  rich  regions. 

The  sparsity  of  data  results  in  large  errors  in  the  divergent  wind  which 
compounds  humidity  errors  through  moisture  convergence.  During  the  initial  hours  of 
model  integration  the  divergent  motion,  diabatic  heating,  and  surface  pressure  fields 
undergo  an  adjustment  in  response  to  the  low  equilibration  of  the  humidity  field.  This 
is  the  period  of  so-called  model  spin-up.  Such  initial  deficiencies  in  the  humidity 
analysis  are  not  addressed  by  the  'dynamical'  type  of  initialization  procedures.  In  this 
context  there  is  a  need  for  enhancement  of  the  information  content  of  the  initial 
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moisture  field.  Physical  initialization  is  a  scheme  that  has  been  developed  to  do  this. 

3.2  Description  of  physical  initialization 

The  term  "physical  initialization"  was  first  introduced  by  Krishnamurti  et  al. 
(1984).  The  details  of  the  initialization  procedure  are  presented  in  Kri^amuiti  et  al. 
(1991).  The  initialization  procedure  includes  the  following  computations: 

1)  diagnostic  measures  of  the  surface  fluxes  of  sensible  heat  and  moisture  according  to 
Yanai  et  al.  (1973), 

2)  generation  of  data  sets  of  the  potential  temperature  and  the  moisture  variable  on  top 
of  the  constant  flux  layer  consistent  with  the  prescribed  fluxes  by  using  the  reverse 
similarity  theory, 

3)  an  analysis  of  the  humidity  variable  above  the  constant  flux  layer  consistent  with 
the  imposed  precipitation  rates  by  the  reverse  cumulus  parameterization,  and 

4)  a  matching  of  the  earth's  radiation  budget  as  inferred  from  the  radiation  algorithm 
and  as  measured  by  the  satellites.  This  procedure  is  completed  by  a  Newtonian 

relaxation  of  all  the  basic  variables  of  the  model  following  a  procedure  outlined  in 
Krishnamurti  et  al.  (1988).  TTiis  procedure  is  schematically  illustrated  in  Figure  3.1. 
In  the  proceeding  sections  the  initialization  scheme  following  Krishnamurti  et  al. 
(1991)  is  described  in  detail. 

3.2.1  Yanai  fluxes 

The  surface  fluxes  of  sensible  heat  (F  )  and  moisture  (F  )  following  Yanai  et  al. 

S 

(1973)  may  be  expressed  by  relations, 

t  A  A 

fJ  =  Qi-Q^-LP 


(3.1) 
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F  J  =  LP  -  ^  (3.2) 

where  Qi,  and  ^  denote  the  vertically  integrated  apparent  heat  source,  net  radiative 
heating,  and  apparent  moisture  sink,  respectively,  and  P  the  observed  rainfall  rate. 

The  apparent  heat  source  and  the  apparent  moisture  sink  are  defined  as 

Qj=^  +  V.sv  +  ^sm  (3.3) 

dt  dp 

02  =  .L(5  +  ^  (3.4) 

dt  dp 

where  s  is  the  large  scale  dry  static  energy. 

3.2.2  Reverse  similarity  theory 

The  similarity  theory  calculates  the  surface  fluxes  of  heat  and  moisture  by  using 
the  similarity  exchange  coefficients  with  the  known  model  variables  (momentum,  heat 
and  moisture)  on  the  top  and  bottom  of  the  surface  layer.  The  reverse  similarity  theory 
calculates  the  temperature  and  moisture  variable  on  the  top  of  surface  layer  by  using 
the  same  similarity  exchange  coefficients  with  the  prescribed  surface  fluxes  of  heat  ana 
moisture. 

Given  the  Yanai  fluxes  of  sensible  and  latent  heat  as  an  input  to  the  similarity 
theory,  we  can  solve  for  the  potential  temperature  and  the  moisture  variable  on  top  of 
the  constant  flux  layer.  Different  sets  of  similarity  equations  describe  the  stable  and 
the  unstable  layers. 

Following  Chang  (1978),  the  similarity  fluxes  of  momentum,  heat  and  moisture 


FORECAST 


Figure  3.1  An  outline  of  the  physical  initialization  within  the  assimilation  phase 
of  the  global  model.  The  package  of  reverse  similarity,  OLR 
matching,  reverse  cumulus  parameterization  slides  from  left  (hour-24) 
to  right  (hour  0)  during  this  initialization  period  (from  Krishnamurti  et 
ai,  1994). 
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are  given  by  the  relations 


•  Ui) 


Fh  =  pCpC„(U2  -  Uixej  -  e,) 


Fq  =  pC,(U2  -  UiXqa  -  qOg^ 


where  the  similarity  exchange  coefficients  are  expressed  as  follows: 


For  the  stable  and  neutral  case:  P4.>0 


c.  = 


u 


k2 


**  <U2  -  Ui?”an(Z2/Zi)]2(l  +  4.7RiJ2 


c„  = 


U  0. 


■k2 


(U2  -  u  -  Si)  0.74[ln(Z2/Zi)]2(l44.7Rig)2 


Cq  =  1.7 


For  the  unstable  case:  R;  <  0 


k2  9.4Ri 

(1 - 


(U2  -  Ui)  [ln(Z2/Z,)]2  1+ClR: 


(3.5) 

(3.6) 

(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


S2 


o  -k2 

^  ~  ________________  —  *■  _ — 1  ) 

"  (Ua  -  U,)(02  -  Oi)  0.74[ln(Z2/Z,)]2  1+C|R: 

o 


(3.12) 


Cq  =  1.7  Cy 


(3.13) 


where. 

Cm*  Cjj,  Cq  =  exchange  coefficient  of  momentum,  heat  and  moisture 
U,  =  friction  velocity 
0.  =  scaling  temperature 

Za,  Zi  =  height  of  top  and  bottom  of  surface  layer 
Ua,  U|  =  wind  at  top  and  bottom  of  surface  layer 
62,  Bi  ss  potential  temperature  at  top  and  bottom  of  surface  layer 
qa,  qi  s  moisture  at  top  and  bottom  of  surface  layer 
ss  ground  wemess 
k  =  von  Karman  constant 

Rig  =  bulk  Richardson  number  which  can  be  expressed  as 


e  (Ua  -  Ui)2 


(3.14) 


The  expression  for  C  follows  the  analysis  of  Louis  (1979): 


C  = 


Cq  X  9.4k^ 
[ln(Z2/Zi)]2  [Z,_ 


(3.15) 


where  Cq  is  taken  to  be  7.4  for  momentum  and  5.3  for  heat  and  moisture. 

To  solve  for  62  and  qa  given  the  surface  fluxes  of  heat  and  moisture,  the  reverse 
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similarity  theory  is  invoked  with  the  assumption  that  the  wind  at  the  lowest  model  level 
U2  is  known.  The  only  unknown  in  the  equations  (3.9)  and  (3.10)  is  62. 

For  the  stable  case,  equation  (3.9)  may  be  written  as 

4.72Ri*  +  (9.4-A)Ri^  +  l=0 

where 

_ V  (U^-Ui)^  k^ 

0.74  gAz  U,e.  [ln(Z2/Zi)l2 

The  signs  of  both  roots  of  the  quadratic  equation  (3.16)  are  positive  and  the  root  with 
the  larger  value  is  beyond  the  physical  upper  limit  of  Ri..  With  Ri  computed.  62  is 
dten  obtained  from  equation  (3.14). 

For  the  unstable  case  (U,6,  <  0)  the  modified  62  is  obtained  by  minimizing  an 
objective  function  using  the  Newton-Raphson  method.  The  objective  function  is 
defined  by 


(3.16) 

(3.17) 


U.0. 

F  = - +  C  (3.18) 

(U2-Ui)(e2-0i)  " 

Computation  of  new  62  is  detailed  in  Krishnamurti  et  al.  (1991).  Once  $2  is  computed, 
C  and  Cq  can  then  be  obtained  using  equations  (3.9),  (3.10),  (3.12)  and  (3.13),  and  q2 
is  obtained  from  equation  (3.7). 

3.2.3  Reverse  cumulus  parameterization 

The  vertical  distribution  of  specific  humidity  is  reanalyzed  such  that  the  rainfall 
implied  by  the  cumulus  parameterization  algorithm  closely  matches  the  prescribed 
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"ot^rved”  rainfall  rates.  The  procedure  for  the  construction  of  a  reverse  cumulus 
parameterization  algorithm  was  discussed  in  Krishnamurti  et  al.  (1988).  This  was 
designed  for  the  modified  Kuo  scheme. 

In  the  reverse  Kuo  (modified),  the  specific  humidity  q  is  modified  as 


Qm 


(3.19) 


where  qm  is  the  modified  specific  humidity,  q  is  the  specific  humidity  prior  to 
modification,  Rq  is  the  observed  rainfall  rate,  q  is  the  vertically  averaged  specific 
humidity  expressed  by 


(3.20) 


The  remaining  symbols  are  explained  in  Chapter  n  (section  2.2). 

The  convergence  of  the  moisture  which  has  been  modified  to  q^  matches  the 
observed  rainfall  Rq.  Thus 

Pt 

-if  m^dp  =  Ro  (3.21) 

®  Jpb  ap 


0 

The  adjustment  procedure  is  such  that  the  total  precipitable  water  in  the  column 
remains  invariant.  That  is. 


Pt  Pt 

-If  qmdp  =  -M  qdp  (3.22) 

^-'Pb  ^■'Pb 

To  obtain  the  analyzed  rainfall,  a  mix  of  raingauge  and  OLR  based  regression 
methods  is  used.  This  method  is  described  in  the  papers,  Krishnarauiti  et  al.  (1983, 
1991,  1993). 

The  reverse  cumulus  parameterization  is  only  applied  to  the  tropical  latitudes 
between  30®S  and  30®N. 

3.2.4  Initialization  of  the  earth's  radiation  budget 

This  initialization  entails  matching  the  OLR  derived  from  radiation  algorithm  of 
the  model  with  the  satellite  observed  OLR.  By  introducing  a  structure  function  (with 
parameter  c)  for  the  moisture  variable  above  500  mb,  the  local  difference  between  two 
estimates  of  OLR  can  be  minimized.  This  minimization  determines  an  optimal  value 
of  €.  This  procedure  tends  to  improve  the  high  and  the  middle  clouds  and  the  planetary 
albedo,  thereby  improving  the  overall  radiation  budget. 

The  minimization  algorithm  utilizes  an  iterative  approach.  For  iteration  i  we 

write 

=  q^*(l  +  £^)  (3.23) 

where  for  f  =  0,  qO  =  q  (analysis). 

The  residue  at  the  end  of  iteration  l-\  is  defined  by 

R'=olr,„-olr' 


(3.24) 
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where  the  subscript  M  denotes  the  model  based  value  and  SAT  denotes  the  satellite 
based  value. 

The  structure  parameter  at  the  iteration  is  computed  by  using  the  bisection 
method.  This  can  be  written  by 


100 


(3.25) 


where  F*’  is  the  convergence 


lauiui  afc  UI&  iKaauuu 


(3.26) 

2 

Here,  an  initial  value  of  the  convergence  factor,  F®,  is  set  to  1. 


3.2.5  Newtonian  relaxation 

This  initialization  is  carried  out  using  a  Newtonian  relaxation  (Hoke  and  Anthes, 
1976)  of  the  basic  model  variables  during  a  preintegration  phase. 

The  equation  for  the  Newtonian  relaxation  takes  the  form, 

^  =  F(A,t)  +  N(A,t)(A®  -  A)  (3.27) 

at 

where  N  denotes  the  relaxation  coefficient.  A®  represents  a  future  value  of  a  variable  A. 
F(A,t)  is  a  forcing  term  of  the  equations  for  variable  A. 

The  integrations  are  carried  out  in  two  steps.  First,  the  tendencies  for  the  normal 
forcing  terms  F  are  computed.  The  next  step  applies  the  Newtonian  term.  This  is 
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expiessed  in  finite  difference  form: 

A(t-f At)  -  A  (t+At)  ^  [AO(t+At)  -  A(t+At)l  (3.28) 

2At 

iRdiere  A*(t+At)  denotes  a  predicted  value  of  A  at  time  (t+At)  prior  to  the  Newtonian 
relaxation.  The  Newtonian  coefficients  are  time  invariant  during  the  integration.  The 
values  used  in  this  study  are  those  used  in  Krishnamurti  et  al.  (1991):  1  x  10^  sec"'  for 
both  vorticity  and  surface  pressure  and  5  x  10“^  sec“*  for  divergence.  The  relaxation 
coefficient  for  humidity  field  is  taken  to  be  5  x  10“^  sec"'  for  the  surface  layer.  No 
relaxation  of  the  humidity  field  is  carried  out  above  the  surface  layer. 

3.3  Some  previous  results  from  physical  initialization 

A  large  number  of  studies  examining  the  impact  of  physical  initialization  have 
been  carried  out  Most  of  these  studies  address  the  impact  on  rainfall  analysis  and 
forecast 

Krishnamurti  et  al.  (1991)  have  illustrated  the  impact  of  physical  initialization  on 
medium  range  forecasts  over  the  tropics.  They  showed  that  rainfall  amounts  and 
locations  of  heavy  rainfall  in  4-5  day  forecasts  from  physically  initialized  initial  fields 
exhibited  close  agreement  with  observed  rainfall.  Mathur  et  al.  (1992)  made  an  effort 
to  study  the  impact  of  physical  initialization  with  NMC  model.  They  demonstrated  a 
strong  positive  impact  on  tropical  prediction  from  the  inclusion  of  a  reverse  cumulus 
parameterization  within  the  six  hourly  assimilation  cycle.  The  improved  initial  state  in 
the  forecast  model  led  to  the  enhancement  of  the  forecast  rain  and  circulation. 
Krishnamurti  et  al.  (1993)  have  shown  the  improvement  in  nowcasting  and  short-range 
forecast  skill  of  rainfall  prediction  by  implementing  the  rainfall  retrieval  algorithm 
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using  OLR,  SSMA,  and  raingauge  measurements  within  the  physical  initialization. 
Kiishnamurti  et  al.  (1994)  explored  the  impact  of  physical  initialization  on  tropical 
rainfall  assimilation  and  forecasts.  They  found  that  the  mean  nowcasting  skill  of 
rainfall  exceeded  the  current  operational  skill. 


CHAPTER  IV 


NUMERICAL  EXPERIMENTS  AND  RESULTS 

4.1  Initial  data 

In  this  study  two  different  types  of  data  have  been  used  as  an  initial  state  for  the 
cloud  forecasts  beginning  at  12  UTC  20  October  1991.  The  first  initial  state  is  taken 
from  ECMWF  nib  analyses.  We  refer  to  this  as  the  normal  or  traditional  initial 
conditions.  The  second  initial  state  consists  of  data  physically  initialized  by  the 
FSUGSM. 

4.1.1  The  normal  initial  data 

The  ECMWF  nib  analyses  provide  the  meteorological  fields  of  height,  zonal  and 
meridional  wind  components,  temperature,  and  relative  humidity  on  12  mandatory 
pressure  levels  (1000,  850,  700,  500,  400,  300,  250,  200,  150,  100,  70,  and  50  mb). 
These  analyses  are  produced  by  the  4-dimensional  data  assimilation  system  (an 
intermittent  system)  which  uses  a  multivariate  optimum  interpolation  analysis,  a 
non-linear  normal  mode  initialization  and  a  high-resolution  forecast  that  produces  the 
first  estimate  for  the  subsequent  analysis.  The  analysis  is  divided  into  two  parts:  one 
for  the  analysis  of  humidity  and  a  second  for  the  simultaneous  analysis  of  geopotential, 
wind,  and  surface  pressure.  A  complete  description  of  the  ECMWF  data  assimilation 
system  is  provided  by  Bengtsson  et  al.  (1982). 
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4.1.2  Physically  initialized  data 

Physical  initialization  was  carried  out  during  a  pre-integration  phase  between  day 
-1  (12  UTC  19  October,  1991)  and  day  0  (12  UTC  20  October,  1991)  using  the  T106 
FSUGSM.  During  the  assimilation  phase,  ingestion  of  the  observed  rainrates  and  OLR 
matching  are  performed  as  shown  in  Figure  3.1.  The  observed  rainrates  are  obtained 
from  an  algorithm  (Krishnamurti  et  al.,  1993)  that  combines  measurements  of 
satellite-based  OLR,  SSM/I,  and  raingauge  data. 

A  modification  of  the  initial  state  variables  via  incorporation  of  analyzed  tropical 
rainrates  is  accomplished  during  the  assimilation  phase  of  the  model  forecast.  The 
physical  initialization  process  produces,  in  a  diagnostic  sense,  a  thermodynamic 
consistency  between  the  humidity  variable,  the  surface  fluxes,  rainfall  distributions, 
diabatic  heating,  and  the  clouds. 

4.2  Results  from  global  forecast  experiments 

Sensitivi^  of  cloud  predictive  skill  was  tested  in  the  FSUGSM  as  a  function  of 
both  initial  conditions  and  implicit  versus  explicit  cloud  schemes.  Four  experiments 
were  designed  as  follows; 

1)  the  experiment  conducted  with  the  implicit  cloud  scheme  using  the  normal  initial 
data  (’contror), 

2)  the  experiment  with  the  same  model  as  in  'control'  but  with  the  physically 
initialized  data  ('EXPl'), 

3)  the  experiment  conducted  with  the  explicit  cloud  scheme  using  the  normal  initial 
data  ('EXP2'),  and 

4)  the  experiment  with  the  same  model  as  in  'EXP2'  but  with  the  physically 
initialized  data  ('EXP3'). 
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A  six  day  forecast  was  made  in  each  of  these  four  experiments. 

4.2.1  Initial  clouds 

4.2.1. 1  Zonal  mean  of  cloudiness 

In  Hgure  4.1  the  zonal-mean  total  cloudiness  for  two  initial  conditions  is 
compared  with  the  surface-based  observations  analyzed  over  land  by  Warren  et  al. 

t 

(1986).  The  ground  data  were  compiled  from  an  eleven  year  data  set  extending  1971  to 
1981.  The  observed  values  (Figure  4.1)  show  a  maximum  cloud  cover  near  the  equator 
and  minimum  in  the  subtropics  of  both  hemispheres.  This  is  a  typical  feature  of  ITCZ’s 
seasonal  position.  Other  peaks  of  maximum  cloud  cover  were  found  near  60  and  60'’N 
which  are  a  result  of  middle  latitude  system  such  as  cold  fronts.  The  initial  cloudiness 
based  on  the  physical  initialization  shows  better  agreement  in  both  the  shape  and 

k) cation  of  the  extreme  values  with  the  observed  data  than  that  of  the  normal  initial 
data.  In  particular,  in  the  tropics  it  shows  a  remarkable  improvement  in  the  amount 
compared  with  the  normal  initial  data.  Obviously  it  is  the  effect  of  physical 
initialization.  Cloud  amounts  of  the  physically  initialized  data  are  consistently  higher 
than  those  of  the  normal  initial  data  in  both  hemispheres. 

4.2. 1.2  Geographical  distribution  of  cloudiness 

Figure  4.2  shows  the  initial  high,  middle,  low  and  total  cloud  distributions  with 
the  normal  data  Oeft)  and  with  the  physically  initialized  data  (right).  A  comparison  of 
the  distribution  of  cloudiness  in  the  two  initial  conditions  shows  that  the  physically 
initialized  data  is  more  representative  in  terms  of 

l)  the  total  cloud  amount, 

2)  the  deep  convective  cloud  associated  with  the  ITCZ  and  Southern  Pacific 
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Convergence  Zone  (SPCZ)  over  active  convective  areas  such  as  the  western 
Pacific  and  northeast  of  Australia,  and  convective  cloud  over  South  America. 

3)  the  high  and  low  cloud  amount  in  the  middle  latitudes  (their  location  matches 
with  the  normal  initial  data), 

4)  the  location  and  amount  of  high  cloud  both  in  the  tropics  and  middle  latitudes 
consistent  with  the  improvement  in  convective  cloud. 


4.2J2  Comparison  of  control  and  EXPl 

4.2.2. 1  Global-mean  cloudiness 


Figure  4.3  shows  the  global-mean  cloud  cover  for  high,  middle,  low  and  total 
clouds  as  a  function  of  forecast  time  for  the  control  run  and  EXPl.  Figure  4.3  (top) 
illustrates  the  fact  that  the  cloud  cover  stabilizes  slowly.  This  suggests  a  substantial 
spin-up  time  for  model  cloud  cover.  This  is  further  illustrated  by  the  following 
tabulation  of  the  spin-up  time  for  high,  middle,  low,  and  total  clouds  derived  from 
Figure  4,3. 


Table  4.1  Percent  cloud  cover  and  spin-up  time  for  control  and  EXPl. 


Initial  Cloud  cover  Spin-up  time 

cloud  cover  after  stabilization  (days) 

(%)  (%) 


Control 


High 

11 

25 

3.5 

Middle 

8 

16 

2.5 

Low 

16 

24 

2.0 

Total 

28 

42 

2.5 

High 

24 

24 

0 

Middle 

17 

17 

0 

Low 

17 

25 

1.5 

Total 

40 

45 

1.5 

EXPl 


cloud  cover  {%) 


zonal  meaui  total  cloud  cover 


Figure  4.1  Comparison  of  the  zonal-mean  total  cloudiness  for  the  normal  initial 
data  (dashed  line)  and  the  physically  initialized  data  (solid  line)  with  a 
surface-based  climatology  (bold  line). 
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Figure  4.2  Global  distribution  of  the  initial  high,  middle,  low  and  total  clouds 
with  the  normal  initial  data  (left)  and  with  the  physically  initialized 
data  (right).  Contour  interval  is  20%.  Values  between  20%  and  80% 
are  dotted;  values  above  80%  are  shaded. 
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Table  4.1  shows  that  the  spin-up  time  for  the  high  clouds  is  the  highest  followed 
by  the  middle  clouds.  Bgure  4.3  (bottom)  shows  that  in  'EXPl'  high  and  middle  clouds 
remain  stable  throughout  the  integration  while  the  low  clouds  tend  towards 
stabilization.  The  spin-up  time  in  EXPl  is  dramatically  less  than  in  control  run.  Also 
note  that  the  total  cloud  cover  in  EXPl  (45%)  is  consistently  higher  than  in  the  control 
run  (42%).  The  former  is  closer  to  the  observational  estimate  (about  50%)  than  the 
latter. 

4.2.2.2  Geographical  distribution  of  cloudiness 

Figures  4.4(a)  through  4.4(f)  illustrate  the  global  cloud  cover  distribution  for  the 
control  run  Geft)  and  EXPl  (right).  Discernible  features  such  as  the  climatological 
cloud  patterns  associated  with  the  ITCZ  and  SPCZ  are  the  criteria  for  comparison 
between  these  experiments.  It  is  seen  from  the  figures  that  EXPl  shows  a  more  robust 
SP<22  over  north  of  Australia  and  Southeast  Asia  and  ITCZ  around  10®N  over  all  6 
days  of  the  forecast  Further,  in  EXPl  the  minimal  cloud  cover  over  Middle  East, 
North  Africa  and  Central  Asia  is  in  consonance  with  climatological  divergent 
circulations  in  the  tropics  (Krishnamurti  et  aL,  1983).  Lastly,  the  bimodality  in  the 
cloud  distribution  over  the  tropical  western  and  central  Pacific  (Ogura  and  Cho,  1973) 
is  much  more  evident  in  EXPl  than  the  control. 

4.2.3  Comparison  of  EXP2  and  EXP3 

4.2.3. 1  Global-mean  cloudiness 

Figure  4.5  shows  the  global-mean  cloud  cover  similar  to  Figure  4.3  for  EXP2  and 
EXP3.  As  before,  a  tabulation  of  the  spin-up  time  derived  from  Figure  4.5  is  presented 
below. 
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Table  4.2  Percent  cloud  cover  and  spin-up  time  for  EXP2  and  EXP3. 


Initial  Cloud  cover  Spin-up  time 

cloud  cover  after  stabilization  (days) 

(^)  (%) 


EXP2 


High 

11 

16 

0.5 

Middle 

8 

12 

2.0 

Low 

16 

30 

3.0 

Total 

28 

39 

2.5 

High 

24 

13 

2.0  (decreased) 

Middle 

17 

14 

1.0  (decreased) 

Low 

17 

30 

2.5 

Total 

40 

41 

2.5 

The  response  of  the  explicit  scheme  to  physically  initialized  data  seems  very 
uxegular.  There  is  an  initial  decay  of  the  high  and  middle  clouds,  and  an  initial  growth 
of  low  clouds  before  stabilization.  This  requires  further  investigation  of  the  explicit 
scheme. 


4.2.3.2  Geographical  distribution  of  cloudiness 

Figure  4.6  illustrates  the  high,  middle,  low  and  total  cloud  distributions  predicted 
by  the  FSUGSM  explicit  scheme  with  normal  initial  data  (left)  and  the  physically 
initialized  data  (right). 

Overall,  high  and  middle  clouds  in  both  experiments,  EXP2  and  EXP3,  tend  to 
decrease  over  the  tropics  during  the  forecast.  E>eep  convective  clouds  are  not  found 
over  the  tropics  in  EXP2,  while  in  EXP3  only  a  few  deep  convective  clouds  are  found 
in  the  western  Pacific.  An  undesirable  feature  in  both  experiments  is  an  increase  in 
cloudiness,  particularly  low  clouds,  during  the  integration.  Therefore,  total  cloud 
amounts  are  mostly  composed  of  the  low  clouds. 
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Figure  4.4(a)  Global  distribution  of  the  high,  middle,  low  and  total  clouds  for  day  1 
forecast  with  the  control  run  Oeft)  and  EXPl  (right).  Contour  interval 
is  20%.  Values  between  20%  and  80%  are  dotted;  values  above  80% 
are  shaded. 
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Figure  4.4(c)  As  Figure  4.4  (a)  but  for  day  3  forecast. 


Figure  4.4(d)  As  Figure  4.4  (a)  but  for  day  4  forecast. 
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Figure  4.6(b)  As  Figure  4.6  (a)  but  for  day  2  forecast. 


Figure  4.6(c)  As  Figure  4.6  (a)  but  for  day  3  forecast 
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Figure  4.6(e)  As  Figure  4,6  (a)  but  for  day  5  forecast. 
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Figure  4.6(f)  As  Figure  4.6  (a)  but  for  day  6  forecast. 
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Beyond  day  3  no  difference  can  be  discerned  in  the  cloud  evolution  between  the 
two  experiments.  In  other  words,  the  impact  of  physical  initialization  on  the  explicit 
cloud  scheme  is  negligible  after  2  days. 


4.2.4  Cloud  predictive  skills 

In  order  to  gauge  the  cloud  forecast  skill  of  the  model,  the  root  mean  square  error 
(RMSE)  is  computed.  We  define  RMSE  as 


RMSEs 


1  Vi 


(Cf  -  C«)2 


where  and  O  denote  the  predicted  and  analyzed  cloud  cover,  respectively,  and  the 
overbar  represents  a  global  average.  In  the  calculation  of  persistence  RMSE  (PRMSE), 
O  is  replaced  by  the  initial  clouds. 

Hgure  4.7  plots  the  ratio  of  RMSE  to  PRMSE  of  the  total  cloud  cover  as  a 
function  of  forecast  time  for  the  control  run,  EXPl,  EXP2  and  EXP3.  This  figure 
clearly  shows  that  EXPl  has  some  predictive  skill  up  to  day  4.  The  control  run 
maintains  the  same  level  of  skill  only  out  to  day  2.5.  Beyond  these  limits  the  RMSE 
ratio  exceeds  1  indicating  that  the  persistence  forecast  is  better.  The  other  experiments 
consistently  have  a  RMSE  ratio  greater  than  1.  This  suggests  that  the  physically 
initialized  initial  state  with  the  implicit  cloud  scheme  has  the  best  predictive  skill  for 
cloud  cover.  All  explicit  experiments  (EXP2,  EXP3)  drift  to  the  climatological  state 
within  day  1  forecast.  The  graph  shows  that  the  explicit  cloud  scheme  has  no 
predictive  skill  for  both  initial  conditions.  This  supports  our  earlier  finding  that  the 
explicit  scheme  has  a  very  irregular  behavior  and  needs  further  investigation. 
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Figure  4.7  The  ratio  of  RMSE  to  PMRSE  of  the  total  cloud  cover  as  function  of 
forecast  time  for  the  control  run  (A),  EXPl(B),  EXP2(C)  and  EXP3 
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Figure  4.8  Scatter  diagram  of  cloud  shortwave  forcing  versus  longwave  forcing 
averaged  for  6  day  forecast  over  the  tropical  Pacific  (120E-120W, 
lOS-lON)  with  the  control  run  (top)  and  EXPl(bottom). 
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Figure  4.9(a)  Distribution  of  the  mean  high,  middle,  low  and  total  clouds  for  6  days 
forecast  over  the  tropical  Pacific  (120E-120W,  lOS-lOI^  with  the 
control  run.  Contour  interval  is  20%.  Values  between  20%  and  80% 
are  shaded;  values  above  80%  are  shaded. 
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linear  fit  is  0.92  for  the  control  run  and  0.91  for  EXPl.  These  underestimate  the  0.96 
correlation  computed  from  ERBE  observations.  The  cloud  forcing  slopes  are  found  to 
be  -2.5  tor  both  experiments.  This  slope  is  steeper  than  that  in  ERBE  observations 
(-1.2).  We  attribute  this  difference  to  the  fact  that  our  values  are  from  a  single  6  day 
forecast  rather  than  from  monthly  mean  observations.  This  discrepancy  also  leads  to  a 
large  difference  with  the  magnitude  of  cloud  forcings  between  the  experiments  and 
observations  i.e.,  360  W/m^  (shortwave),  150  W/m^  (longwave)  for  both  experiments; 
150  W/m2  (shortwave),  120  W/m^  (longwave)  for  observations.  There  are  a  few  subtle 
differences  that  can  be  discerned  from  Figure  4.8  between  the  control  and  EXPl.  In 
Figure  4.8  (top)  note  the  small  conglomeration  of  points  over  relatively  small  negative 
shortwave  forcing  ranging  -144  to  -288  W/m^  and  relatively  large  longwave  forcing 
ranging  90  to  120  W/m2.  This  suggests  a  significant  amount  of  low  clouds  in  the 
control  experiment  forecast  This  is  substantiated  by  Figure  4.9(a).  Similarly,  in 
Figure  4.8  (bottom)  notice  that  there  is  a  similar  kind  of  conglomeration  of  points  as 
before  but  over  relatively  small  negative  shortwave  forcing  ranging  0  to  144  W/m2  and 
relatively  small  longwave  forcing  ranging  0  to  30  W/m^.  This  indicates  that  EXPl  has 
more  high  clouds  in  its  forecast  Figure  4.9(b)  supports  this  contention. 


CHAPTER  V 


CONCLUSIONS 

In  this  study  extensive  work  on  cloud  forecasts  has  been  carried  out.  Two 
distinct  cloud  parameterizations  have  been  examined  in  terms  of  both 
representativeness  and  predictive  skill.  One  of  these  is  the  classical  implicit  method 
and  the  other  an  explicit  method  which  carries  certain  cloud  variables  as  prognostic 
variables.  For  each  method  two  forecasts  were  made:  one  from  physically  initialized 
data  and  the  other  from  initial  conditions  not  subjected  to  physical  initialization.  These 
experiments  provide  the  basis  for  an  evaluation  of  the  implicit  and  explicit  schemes. 
From  these  experiments  we  also  assess  the  impact  of  the  initial  data  on  the  cloud 
forecast  Based  upon  results  from  the  four  experiments  the  following  conclusions  can 
be  drawn: 

1)  The  initial  cloudiness  of  the  physically  initialized  data  set  was  in  close  agreement 
with  the  observed  cloudiness.  In  the  tropics  it  shows  a  remarkable  improvement 
in  the  cloud  cover.  Cloud  amounts  of  the  physically  initialized  data  are 
consistently  higher  than  those  of  the  normal  data  set  in  both  hemispheres.  The 
distribution  of  initial  cloudiness  shows  that  the  physically  initialized  data  is  more 
representative  in  terms  of  a)  the  total  cloud  amount,  b)  the  deep  convective  cloud 
associated  with  ITCZ  and  SPCZ,  over  active  convective  areas  such  as  the  western 
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Pacific,  northeast  of  Australia,  and  over  northern  part  of  South  America,  and  c) 
the  high  and  low  cloud  amount  in  the  middle  latitudes. 

2)  The  physically  initialized  data  set  has  shown  a  major  reduction  in  the  spin-up 
time  of  high,  middle  and  low  clouds  with  the  implicit  cloud  scheme. 

3)  The  response  of  the  explicit  scheme  to  physically  initialized  data  was  very 
irregular.  There  was  an  initial  decay  of  the  high  and  middle  clouds  and  an  initial 
growth  of  low  clouds  before  stabilization.  Overall,  high  and  middle  clouds  in 
both  explicit  experiments  tend  to  decrease  over  the  tropics  during  the  forecast. 
Beyond  day  3  no  difference  is  noticed  in  the  cloud  evolution  between  the 
experiments  with  and  without  the  physical  initialization.  Detailed  study  of  the 
explicit  cloud  scheme  is  an  avenue  for  further  research. 

4)  The  physically  initialized  initial  state  with  the  implicit  cloud  scheme  has  the  best 
predictive  skill  (up  to  day  4).  The  normal  initial  state  with  the  implicit  cloud 
scheme  maintains  the  same  level  of  skill  only  up  to  day  2.5.  All  explicit 
experiments  drift  to  the  climatological  state  within  day  1  of  the  forecast 

5)  There  is  a  strong  negative  correlation  between  the  cloud  shortwave  forcing  and 
longwave  forcing  (i.e.,  -  0.92  for  the  normal  initial  data  and  -  0.91  for  the 
physically  initialized  data).  This  is  in  good  agreement  with  the  ERBE 
observations  of  -  0.96.  Furthermore  it  is  found  that  the  implicit  scheme  without 
the  physical  initialization  tends  to  produce  a  larger  population  of  low  clouds 
(small  negative  shortwave  forcing,  large  longwave  forcing),  while  that  with  the 
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physical  initialization  produces  a  larger  population  of  high  clouds  (small  negative 
shortwave  forcing,  small  longwave  forcing). 
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